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with  a  modulation  frequency  which  is  equal  to  the  fundamental  vibrational  frequency 
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I.  FORWARD 


The  central  theme  of  our  work  over  the  past  three  years  has  been  the  use  of 
electromagnetically  induced  transparency  (EIT)  for  new  types  of  nonlinear  interactions 
and  processes.  The  last  several  years  have  been  unusually  productive.  In  a 
February  1999  Nature  article  which  received  world- wide  attention,  Lene  Hau,  her 
students,  and  the  author  of  this  report,  Steve  Harris,  showed  that,  by  using  EIT  in  an 
ultra-cold  gas  of  Na  atoms,  an  optical  pulse  could  be  caused  to  travel  107  times  slower 
than  the  speed  of  light  in  vacuum.  In  a  June  1991  Physical  Review  Letter ,  Harris  and 
Hau  showed  how  this  slow  group  velocity  could  result  in  nonlinear  coefficients  which  are 
sufficiently  large  that  it  may  soon  be  possible  to  do  nonlinear  optics  at  picowatt  power 
levels.  In  a  May  1999  Optics  Letter,  the  technique  of  nonlinear  optics  at  maximum 
coherence  was  used  to  demonstrate  efficient  gas  phase  frequency  conversion  to  the 
vacuum  ultraviolet.  In  particular,  a  233-nm  laser  beam  was  converted  to  186  nm  at 
approximately  unity  photon  efficiency.  In  follow-up  work,  this  same  technique  was 
extended  to  a  system  where  all  frequencies  are  resonant.  Here,  a  conversion  efficiency  of 
about  30%  was  obtained  at  laser  power  densities  of  about  105  W/cm2  and  at  atom 
densities  of  only  1014  atoms/cm3.  In  another  area  of  substantial  progress,  Harris  and 
Sokolov  have  suggested  a  new  technique  for  producing  subfemtosecond  and  even  single¬ 
cycle  optical  pulses.  The  idea  is  to  use  two  laser  beams  whose  frequency  difference  is 
slightly  offset  from  a  molecular  transition.  This  produces  a  single,  molecular, 
superposition  state  which,  in  the  manner  of  EIT,  may  be  driven  at  maximum  coherence. 
In  recent  experimental  work,  we  have  shown  the  collinear  generation  of  equidistant 
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sidebands  covering  about  50,000  cm'1  of  spectral  bandwidth  and  ranging  from  2.94  p  in 
the  infrared  to  195  nm  in  the  ultraviolet.  This  work  is  the  first  demonstration  of  a 
nonlinear  optical  process  at  maximum  coherence  in  a  molecular  system.  Because  these 
sidbands  comprise  an  FM  signal,  group  velocity  dispersion  may  be  used  to  temporally 
compress  this  waveform  to  a  train  of  subfemtosecond  pulses. 


II.  STATEMENT  OF  PROBLEM  STUDIED 


We  study  electromagnetically  induced  transparency  (EIT),  both  in  its  own  right  and  as  a 
vehicle  to  make  new  types  of  nonlinear  optical  devices.  Examples  of  such  devices 
include  subfemtosecond  and  single-cycle  pulse  generators  and  nonlinear  optical 
converters  which  may  function  at  the  single-photon  level. 
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III.  SUMMARY  OF  MOST  IMPORTANT  RESULTS 


To  show  the  overall  evolution  of  this  work  the  Summary  of  Important  Results 
includes  the  main  contributions  of  our  earlier  contracts  in  this  area. 

(1)  Our  work  in  this  field  began  with  the  paper:  S.  E.  Harris,  J.  E.  Field,  and 
A.  Imamoglu,  "Nonlinear  Optical  Processes  Using  Electromagnetically  Induced 
Transparency,"  Phys.  Rev.  Lett.  64,  1 107-1 1 10  (March  1990).  This  work  also  introduced 
the  expression  "electromagnetically  induced  transparency." 

(2)  The  first  observation  of  EIT  was  reported  in:  K.-J.  Boiler,  A.  Imamoglu,  and 
S.  E.  Harris,  "Observation  of  Electromagnetically  Induced  Transparency,"  Phys.  Rev. 
Lett.  66,  2593-2596  (May  1991). 

(3)  The  concept  of  matched  pulses  and  the  normal  modes  of  EIT  were  introduced 
in:  S.  E.  Harris,  "Electromagnetically  Induced  Transparency  with  Matched  Pulses," 
Phys.  Rev.  Lett.  70,  552-555  (February  1993)  and  S.  E.  Harris,  "Normal  Modes  for 
Electromagnetically  Induced  Transparency,"  Phys.  Rev.  Lett.  72,  52-55  (January  1994). 

(4)  In  March  of  1995  we  showed  how  EIT  would  cause  well-shaped  pulses  to 
propagate  at  a  group  velocity  as  slow  as  c/165:  A.  Kasapi,  M.  Jain,  G.  Y.  Yin,  and 
S.  E.  Harris,  "Electromagnetically  Induced  Transparency:  Propagation  Dynamics,"  Phys. 
Rev.  Lett.  74,  1447-2450  (March  1995). 

(5)  The  observation  that  EIT  would  also  eliminate  self-focusing  and  defocusing 
was  described  in:  M.  Jain,  A.  J.  Merriam,  A.  Kasapi,  G.  Y.  Yin,  and  S.  E.  Harris, 
"Elimination  of  Optical  Self-Focusing  by  Population  Trapping,"  Phys.  Rev.  Lett.  75, 
4385-4388  (December  1995). 
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(6)  The  idea  of  nonlinear  optics  at  maximum  coherence  was  introduced  and 
experimentally  demonstrated  in  the  paper:  M.  Jain,  H.  Xia,  G.  Y.  Yin,  A.  J.  Merriam, 
and  S.  E.  Harris,  "Efficient  Nonlinear  Frequency  Conversion  with  Maximal  Atomic 
Coherence,"  Phys.  Rev.  Lett.  77,  4326-4329  (November  1996). 

(7)  The  suggestion  for  subfemtosecond  pulse  generation  by  molecular  modulation 
was  made  in:  S.  E.  Harris  and  A.  V.  Sokolov,  "Subfemtosecond  Pulse  Generation  by 
Molecular  Modulation,"  Phys.  Rev.  Lett.  81,  2894-2897  (October  1998). 

(8)  A  suggestion  for  nonlinear  optics  at  near  the  single-photon  level  was  described 
in:  S.  E.  Harris  and  Y.  Yamamoto,  "Photon  Switching  by  Quantum  Interference,"  Phys. 
Rev.  Lett.  81,  3611-3614  (October  1998). 

(9)  The  report  of  slow  light  was  made  in:  L.  V.  Hau,  S.  E.  Harris,  Z.  Dutton,  and 
C.  H.  Behroozi,  "Light  Speed  Reduction  to  17  Meters  Per  Second  in  an  Ultracold  Atomic 
Gas,"  Nature  397,  594-598  (February  1999). 

(10)  The  use  of  EIT  to  efficiently  generate  vacuum  ultraviolet  radiation  was 
described  in:  A.  J.  Merriam,  S.  J.  Sharpe,  H.  Xia,  D.  Manuszak,  G.  Y.  Yin,  and 
S.  E.  Harris,  "Efficient  Gas-Phase  Generation  of  Coherent  Vacuum  Ultraviolet 
Radiation,"  Opt.  Lett.  24,  625-627  (May  1999). 

(11)  The  idea  that  EIT  could  be  obtained  in  a  medium  with  hyperfine  structure 
was  experimentally  shown  in:  H.  Xia,  A.  J.  Merriam,  S.  J.  Sharpe,  G.  Y.  Yin,  and 
S.  E.  Harris,  "Electromagnetically  Induced  Transparency  with  Spectator  Momenta," 
Phys.  Rev.  A  59,  R3190-R3193  (May  1999). 
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(12)  A  description  of  nonlinear  optics  with  slow  light  was  given  in:  S.  E.  Harris 
and  L.  V.  Hau,  "Nonlinear  Optics  at  Low  Light  Levels,"  Phys.  Rev.  Lett.  82,  4611-4614 
(June  1999). 

(13)  A  general  extension  of  EIT  to  greatly  detuned  dispersive  media  was  given  in: 
D.  D.  Yavuz,  A.  V.  Sokolov,  and  S.  E.  Harris,  "Eigenvectors  of  a  Raman  Medium,"  Phys. 
Rev.  Lett.  84, 75-78  (January  2000). 

(14)  We  have  experimentally  demonstrated  the  use  of  molecular  deuterium  in  two 
driving  lasers  to  demonstrate  collinearly  generation  of  mutually-coherent,  equidistant 
sidebands  covering  50,000  cm'1  of  spectral  bandwidth:  A.  V.  Sokolov,  D.  R.  Walker, 
D.  D.  Yavuz,  G.  Y.  Yin,  and  S.  E.  Harris,  "Raman  Generation  by  Phased  and  Antiphased 
Molecular  States,"  Phys.  Rev.  Lett.  85,  562-565  (July  2000). 

(15)  EIT  has  been  used  to  produce  efficient  nonlinear  frequency  in  an  all-resonant 
double-A  system.  At  power  densities  which  are  characteristic  of  cw  lasers,  we  obtain 

30%  conversion  efficiency  from  233  nm  to  the  vacuum  ultraviolet  at  186  nm. 
[A.  J.  Merriam,  S.  J.  Sharpe,  M.  Shverdin,  D.  Manuszak,  G.  Y.  Yin,  and  S.  E.  Harris, 
"Efficient  Nonlinear  Frequency  Conversion  in  an  All-Resonant  Double-L  System,"  Phys. 
Rev.  Lett.  84,  5308-5311  (June  2000).] 

(16)  A  special  feature  of  slowly-moving  light  is  the  associated  spatial 
compression.  If  an  optical  pulse  is  slowed  down  in  the  ratio  of  c/Vg,  it  is  also  spatially 
compressed  by  this  same  factor.  This  compression,  in  turn,  leads  to  extremely  large 
longitudinal  gradient  forces.  The  forces,  in  turn,  result  in  new  types  of  optical  scattering 
as  well  as  ballistic  atomic  motion  and  atom  surfing.  [S.  E.  Harris,  "Ondermotive  Forces 
with  Slow  Light,"  Phys.  Rev.  Lett.  85, 4032-4035  (November  2000).] 
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(17)  In  work  which  will  soon  be  published  we  have  shown  that  the  sidebands 
which  are  obtained  from  our  new  type  of  Raman  generator  may  be  phased  so  as  to 
produce  amplitude-  or  frequency-modulated  light.  The  modulation  frequency  is  equal  to 
the  fundamental  vibrational  frequency  of  molecular  deuterium  (2994  cm'1  ~  90  THz). 
[A.  V.  Sokolov,  D.  D.  Yavuz,  D.  R.  Walker,  G.  Y.  Yin,  and  S.  E.  Harris,  "Light 
Modulation  at  Molecular  Frequencies,"  Phys.  Rev.  A  (to  be  published).] 

We  attach  abstracts  of  all  papers  which  were  published  during  this  report  period 
as  Appendix  A.  Appendix  B  includes  copies  of  four  of  the  most  significant  publications. 
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IV.  RELATION  TO  OTHER  CONTRACTS 


We  note  that  the  work  reported  here  was  jointly  supported  by  the  U.S.  Air  Force 
Office  of  Scientific  Research,  the  U.S.  Office  of  Naval  Research,  and  the  U.S.  Army 
Research  Office  OSD  Multidisciplinary  University  Research  Initiative  Program  (MURI). 
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We  suggest  a  technique  for  producing  subfemtosecond  pulses  of  radiation.  The  technique  is  based 
on  using  ciectromagnetically  induced  transparency  to  produce  a  strongly  driven  molecular  coherence. 
This  coherence  results  in  a  Raman  spectrum  with  Bessel  function  amplitudes  and  phases  corresponding 
to  a  frequency  modulated  signal  thereby  allowing  compression  by  the  group  velocity  dispersion  of  the 
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We  describe  a  tour-state  atomic  system  that  absorbs  two  photons,  but  does  not  absorb  one  photon. 
As  a  switch,  in  the  ideal  limit,  it  operates  at  an  energy  cost  of  one  photon  per  switching  event. 
[S003 1-9007(98)0742 1-3] 

PACS  numbers:  32.80.— t.  42.50. Dv.  42.50. Gy.  42.50. Hz 


letters  to  nature 


Light  speed  reduction 
to  17  metres  per  second 
in  an  ultracold  atomic  gas 

Lene  Vestergaard  Hau't,  S.  E.  Harris*,  Zachary  Dutton** 

&  Cyrus  H.  Behroozi*§ 

*  Rowland  Institute  for  Science,  100  Edwin  H.  Land  Boulevard,  Cambridge, 
Massachusetts  02142,  USA 

*  Department  of  Physics.  $  Division  of  Engineering  and  Applied  Sciences. 
Harvard  University.  Cambridge ,  Massachusetts  02138.  USA 

X  Edward  L.  Ginzton  Laboratory,  Stanford  University ;  Stanford ,  California  94305, 
USA 


Techniques  that  use  quantum  interference  effects  are  being 
actively  investigated  to  manipulate  the  optical  properties  of 
quantum  systems1.  One  such  example  is  electromagnet ically 
induced  transparency,  a  quantum  effect  that  permits  the  propaga¬ 
tion  of  light  pulses  through  an  otherwise  opaque  medium2  5.  Here 
we  report  an  experimental  demonstration  of  electromagnetically 
induced  transparency  in  an  ultracold  gas  of  sodium  atoms,  in 
which  the  optical  pulses  propagate  at  twenty  million  times  slower 
than  the  speed  of  light  in  a  vacuum.  The  gas  is  cooled  to 
nanokelvin  temperatures  by  laser  and  evaporative  cooling**10. 
The  quantum  interference  controlling  the  optical  properties  of 
the  medium  is  set  up  by  a  ‘coupling  laser  beam  propagating  at  a 
right  angle  to  the  pulsed  ‘probe1  beam.  At  nanokelvin  tempera¬ 
tures,  the  variation  of  refractive  index  with  probe  frequency  can 
be  made  very  steep.  In  conjunction  with  the  high  atomic  density, 
this  results  in  the  exceptionally  low  light  speeds  observed.  By 
cooling  the  cloud  below  the  transition  temperature  for  Bose- 
Einstein  condensation1 1"13  (causing  a  macroscopic  population  of 
alkali  atoms  in  the  quantum  ground  state  of  the  confining 
potential),  we  observe  even  lower  pulse  propagation  velocities 
(17ms1)  owing  to  the  increased  atom  density.  We  report  an 
inferred  nonlinear  refractive  index  of  0.18cm:  W"1  and  find  that 
the  system  shows  exceptionally  large  optical  nonlinearities,  which 
are  of  potential  fundamental  and  technological  interest  for  quantum 
optics. 
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Subfemtosecond  pulse  generation  by  rotational 
molecular  modulation 
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Edward  L.  Ginzton  Laboratory,  Stanford  University,  Stanford,  California  94305 
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We  extend  a  recent  suggestion  for  the  generation  of  subfemtosecond  pulses 

[Fhvs  Rev  Lett.  81,  2894  (1998)1  to  the  rotational  spectrum  of  molecular  hydrogen  (H,).  When  a  rotational 
transition  ia>  -  |6>  is  strongly  driven  (|Po4|  -  0.5)  the  generation  and  phase-slip  lengths  are  of  the  same  order 
anTthe  Raman  spectrum  has  approximately  Bessel  function  sideband  amplitudes.  Numerical  simulation 
predicts  that  this  spectrum  (generated  in  a  14-cm-long  cell  at  1-atm  pressure  of  H,)  wiU  compress  mto  a  tram 
of  pulses  with  94-fs  pulse  separation  and  a  pulse  length  of  0.5  fs.  ©  1999  Optical  Society  of  America 
OCIS  codes :  320.0320,  320.5520,  190.3100,  030.1640. 
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Efficient  gas-phase  generation  of  coherent  vacuum 

ultraviolet  radiation 
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Edward  L  Ginzton  Laboratory,  Stanford  University,  Stanford.  California  94305 

Received  February  1,  1999 

We  report  the  demonstration  of  a  pulsed  atomic  lead  (Pb)  vapor-based  vacuum  ultraviolet  frequency  converter 
from  233  to  186  nm  with  unity  photon-conversion  efficiency.  This  conversion  is  attained  without  phase 
matching.  ©  1999  Optical  Society  of  America 
OCIS  codes:  190.2620,  190.4410,  270.1670,  020.1670,  140.7220,  140.4380. 
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Electromagnetically  induced  transparency  with  spectator  momenta 

Hui  Xia,  A.  J.  Merriam,  S.  J.  Sharpe,  G.  Y.  Yin,  and  S.  E.  Harris 
Edward  L  Ginzton  Laboratory ,  Stanford  University,  Stanford,  California  94305 
(Received  18  May  1998) 

We  describe  a  method  and  present  experimental  results  for  obtaining  electromagnetically  induced  transpar¬ 
ency  by  adjusting  the  frequencies  of  two  lasers  so  that  they  coincide  with  the  centers  of  gravity  of  two 
hyperfine  split  transitions.  No  hyperfine  states  are  in  Raman  resonance  and  the  zero  of  the  dipole  moment 
results  from  the  interference  of  the  manifold  of  cooperating  m  states.  [S  1050-2947(99)50805-6] 
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Nonlinear  Optics  at  Low  Light  Levels 

S.E.  Harris 
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Lene  Vestergaard  Hau 

Rowland  Institute  for  Science,  100  Edwin  H.  Land  Boulevard,  Cambridge ,  Massachusetts  02142 
and  Department  of  Physics,  Hansard  University,  Cambridge,  Massachusetts  02138 
(Received  21  December  1998) 

We  show  how  the  combination  of  electromagnetically  induced  transparency  based  nonlinear  optics 
and  cold  atom  technology,  under  conditions  of  ultraslow  light  propagation,  allows  nonlinear  processes 
at  energies  of  a  few  photons  per  atomic  cross  section.  [S003l-9007(99)09290-X] 

PACS  numbers:  32.80.-t.  42.50.Dv.  42.50.Gy,  42.50.Hz 


PHYSICAL  REVIEW  A 


VOLUME  60,  NUMBER  2 


AUGUST  1999 


Subfemtosecond  pulse  generation  with  molecular  coherence  control 
in  stimulated  Raman  scattering 
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A.  V.  Sokolov 

Edward  L  Ginzton  Laboratory,  Stanford  University,  Stanford,  California  94305 
(Received  4  March  1999) 

We  study  the  generation  of  subfemtosecond  pulses  with  the  molecular  coherence  control  in  stimulated 
Raman  scattering.  We  show  analytically  that  the  antiphased  state  temporally  advances  the  higher  frequencies 
with  respect  to  the  lower  frequencies  during  a  beating  cycle.  After  some  propagation  distance,  due  to  the 
dispersion  and  the  difference  between  the  antiphased  and  phased  states  in  advancing  high  or  low  frequencies, 
the  coherence  is  highest  on  the  negative  side  of  the  detuning,  when  the  two-photon  Rabi  frequency  is  about 
equal  to  the  detuning.  This  asymmetry  of  the  coherence  magnitude  with  respect  to  the  negative  and  positive 
sides  of  the  Raman  detuning  is  reflected  in  the  behavior  of  the  tuning  characteristics  of  the  high-order  fre¬ 
quency  components.  When  the  Raman  detuning  is  small,  although  the  process  is  nonadiabatic,  the  subfemto¬ 
second  pulse  generation  may  occur  for  both  negative  and  positive  sides  of  the  detuning. 

[S  1050-2947(99)00908-7] 
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Subfemtosecond  Compression  of  Periodic  Laser  Pulses 
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Edward  L.  Ginzton  Laboratory,  Stanford  University,  Stanford.  California  94305 

(April  23.  1999) 


Abstract 


We  show  that  a  recently  proposed  technique  of  molecular  modulation  can 
be  used  to  compress  mode-locked  laser  pulses  into  the  subfemtosecond  do¬ 
main.  The  method  relies  on  driving  a  molecular  oscillation  with  a  periodic 
pulse  train.  The  coherent  molecular  motion,  in  turn,  modulates  the  laser  fre¬ 
quency  and  allows  pulse  compression  by  group  velocity  dispersion  of  the  same 
molecular  medium. 
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Efficient  Gas-Phase  VUV  Frequency  Up-Conversion 

Andrew  J.  Meniam,  Scott  J.  Sharpe,  Hui  Xia,  Danielle  A.  Manuszak, 

G.  Y.  Yin,  and  Stephen  E.  Harris,  Fellow,  IEEE 


Abstract — We  describe  a  pulsed,  atomic  lead  (Pb)  vapor-based 
vacuum  ultraviolet  (VUV)  frequency  converter  from  233  to  186 
nm  with  near-unity  photon  conversion  efficiency.  This  conver¬ 
sion  is  attained  without  phasematching  and  is  accomplished 
by  using  electromagnetically  induced  transparency  to  drive  a 
Raman  transition  to  near-maximal  coherence  without  loss  or 
beam  biow-up.  Under  these  conditions,  the  linear  and  nonlinear 
polarizations  of  the  generated  186-nm  field  are  of  the  same  order 
and  complete  conversion  from  the  233-nm  field  occurs  within  a 
single  (nonphasematched)  coherence  length. 

Index  Terms — Atomic  physics,  lead,  nonlinearities,  nonlinear 
optics,  optical  frequency  conversion,  optical  parametric  ampli¬ 
fiers,  optical  phasematching. 
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Eigenvectors  of  a  Raman  Medium 

D.  D.  Yavuz.  A.  V.  Sokolov,  and  S.  E.  Harris 

Edward  L.  Ginzton  Laboratory,  Stanford  University,  Stanford,  California  94305 
(Received  18  August  1999) 

We  show  the  existence  of  discrete  sets  of  Raman  sidebands  which  self-consistently  establish  a  Raman 
coherence  and  propagate  without  change  in  amplitude  and  relative  phase.  Equivalently,  there  exist  peri¬ 
odic  femtosecond-time-scale,  temporal  pulse  shapes  which  propagate  without  change  in  shape. 


PACS  numbers:  42.50.Gy,  32.80.Qk.  42.65.Dr.  42.65.Tg 
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Raman  Generation  by  Phased  and  Antiphased  Molecular  States 

A.  V.  Sokolov.  D.  R.  Walker.  D.  D.  Yavuz.  G.  Y.  Yin.  and  S.  E.  Harris 

Edward  L.  Ginzton  Laboratory.  Stanford  University.  Stanford.  California  94305 
(Received  14  December  1999) 

We  use  molecular  deuterium  and  two  driving  lasers  to  demonstrate  collinear  generation  of  mutually 
coherent  equidistant  sidebands,  covering  50000  cm’1  of  spectral  bandwidth  and  ranging  from  2.94  /xm 
to  195  nm  in  wavelength.  The  essential  idea  is  the  adiabatic  preparation  of  a  single,  highly  coherent 
(\p,,h\  -  0.33)  molecular  eigenstate. 


PACS  numbers:  42. 50. Gy.  32.S0.Qk.  42.60. Fc.  42.65. Dr 
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Efficient  Nonlinear  Frequency  Conversion  in  an  All-Resonant  Double- A  System 

Andrew  J.  Merriam.  S.  J.  Sharpe.  M.  Shverdin,  D.  Manuszak.  G.  Y.  Yin.  and  S.  E.  Hams 

Edward  L.  Ginzton  Laboratory ,  Stanford  University.  Stanford.  California  94305 
(Received  3  February  2000) 

We  demonstrate  efficient,  pulsed,  gas-phase,  nonlinear  frequency  conversion  in  a  quadruply  resonant, 
double-A  system  and.  simultaneously,  verify  theoretical  predictions  of  Rabi-frequency  matching  unique 
to  absorbing  nonlinear  media.  This  system  is  used  to  up-convert  ultraviolet  light  at  233  nm  to  the 
vacuum  ultraviolet  at  186  nm  in  atomic  Pb  vapor  with  small-signal  conversion  efficiencies  exceeding 
30%  and  with  modest  atomic  density-length  (NL)  products  (scale  1014  cm-2)  and  optical  power  densities 
(10-100  kW/cm2). 

PACS  numbers:  42.50.Gy.  42.50.Hz.  42.65. Kv.  42.79. Nv 


Pondermotive  Forces  with  Slow  Light 


S.  E.  Harris 

Edward  L.  Ginzton  Laboratory ,  Stanford  University,  Stanford .  California  94305 

(7/14/00) 

Abstract 

This  work  describes  atomic  processes  which  result  from  the  greatly  enhanced 
longitudinal  gradient  force  which  is  inherent  to  the  propagation  of  slow  light. 
These  processes  are:  (1)  ballistic  atom  motion  and  atom  surfing  and  (2)  a 
type  of  local  pondermotive  nonlinearity  or  scattering  which  results  from  free- 

particle  sinusoidal  motion  and  the  density  variation  caused  by  this  motion. 
PACS  numbers:  32.80.~t.  32.80.Lg,  32.80.Qk,  42.50.Vk 
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Collinear  Light  Scattering  Using 
Electromagnetically  Induced 
Transparency 

S.  E.  Harris.  A.  V.  Sokolov.  D.  R.  Walker. 

D.  D.  Yavuz.  and  G.  Y.  Yin 


Edward  L.  Ginzton  Laboratory, 

Stanford  University,  Stanford.  California  94305 


Abstract.  The  paper  describes  two  types  of  nonlinear  optical  processes  which  are 
based  on  electromagnetically  induced  transparency.  These  are:  (1)  Collinear  generation 
of  FM-like  Raman  sidebands  and  (2)  a  type  of  pondermotive  light  scattering  which  is 
inherent  to  the  interaction  of  slow  light  with  cold  atoms.  Connections  to  other  areas 
of  EIT-based  nonlinear  optics  are  also  described. 


Light  Modulation  at  Molecular  Frequencies 


A.  V.  Sokolov.  D.  D.  Yavuz.  D.  R.  Walker.  G.  Y.  Yin.  and  S.  E.  Harris 

Edward  L.  Ginzton  Laboratory.  Stanford  University,  Stanford.  California  94305 

(June  2.  2000) 


Abstract 


We  describe  the  generation  of  amplitude  and  frequency  modulated  light 
with  a  modulation  frequency  equal  to  the  fundamental  vibrational  frequency 
of  molecular  deuterium  (2994  cm-1  ~  90  THz).  This  is  a  first  step  toward 
the  synthesis  of  subfemtosecond  pulses  with  prescribed  temporal  shape. 


PACS  numbers:  42.50.Gy,  32.80.Qk,  42.60.Fc,  42.65.Dr,  42.65.Re.  42.65.Sf 


Optical  Frequency  Conversion  by  a  Rotating  Molecular  Waveplate 


A.  V.  Sokolov.  S.J.  Sharpe.  M.  Shverdin.  D.  R.  Walker.  D.  D.  Yavuz.  G.  Y.  \  in.  and  S.  E. 

Harris 

Edward  L.  Ginzton  Laboratory.  Stanford  University,  Stanford.  California  94305 

(September  6,  2000) 


Abstract 


We  demonstrate  efficient  laser  frequency  conversion  in  low-pressure  molec¬ 
ular  deuterium.  We  use  two  laser  fields,  with  opposite  circular  polarizations, 
to  produce  a  strong  excitation  of  a  ro-vibrational  transition,  at  a  frequency  of 
3167  cm-1.  The  coherent  molecular  motion,  in  turn,  modulates  a  third  laser 
field  (also  circularly  polarized)  and  results  in  highly  efficient  single-sideband 
conversion. 


OCIS  numbers:  190.2620.  190.5650,  270.1670.  300.6270 


Refractive  Index  of  a  Strongly  Driven  Raman  Medium 


A.  V.  Sokolov,  D.  D.  Yavuz,  and  S.  E.  Harris 
Edward  L.  Ginzton  Laboratory,  Stanford  University,  Stanford,  California  94.305 

(December  14,  2000) 


Abstract 


We  derive  an  expression  for  the  refractive  index  and  self-focusing  (or  de- 
focusing)  distance  in  a  strongly  driven  Raman  medium,  when  the  coherence 
of  the  molecular  transition  is  on  the  order  of  its  maximum  value  |p«*l  ~  °-5- 
In  this  regime  Raman  generation  and  sideband  propagation  are  intrinsically 
interrelated,  and  the  refractive  index  is  determined  by  the  strong  interaction 

of  the  whole  comb  of  Raman  sidebands  with  the  medium. 

PACS  numbers:  42.65.Re,  32.80.Qk,  42.50.Gy,  42.65.Sf 
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Subfemtosecond  Pulse  Generation  by  Molecular  Modulation 

S.  E.  Harris  and  A.  V.  Sokolov 

Edward  L.  Ginzton  Laboratorx.  Stanford  University,  Stanford ,  California  94305 
(Received  18  May  1998) 

We  suggest  a  technique  for  producing  subfemtosecond  pulses  of  radiation.  The  technique  is  based 
on  using  electromagnetically  induced  transparency  to  produce  a  strongly  driven  molecular  coherence. 
This  coherence  results  in  a  Raman  spectrum  with  Bessel  function  amplitudes  and  phases  corresponding 
to  a  frequency  modulated  signal,  thereby  allowing  compression  by  the  group  velocity  dispersion  of  the 
same  medium.  [S003 1-9007(98)07282-2] 


PACS  numbers:  42.65.Re,  32.80.Qk,  42.50.Gy,  42.65.Sf 

This  Letter  suggests  and  analyzes  a  technique  for  pro¬ 
ducing  subfemtosecond  pulses  of  radiation.  We  will  show 
that  two  laser  beams  whose  frequency  difference  is  slightly 
offset  from  a  molecular  transition  will,  for  an  appropriate 
choice  of  gas  pressure  and  cell  length,  generate  a  spec¬ 
trum  of  Raman  sidebands  whose  Fourier  transform  is  a 
periodic  train  of  subfemtosecond  pulses.  The  essence  of 
the  technique  is  the  concurrent  generation  of  a  frequency 
modulated  (FM)  waveform  and  the  use  of  group  velocity 
dispersion  to  temporally  compress  this  waveform.  By  nu¬ 
merically  modeling  this  process  in  molecular  deuterium 
{Df),  we  calculate  a  generated  train  of  pulses  with  a  pulse 
spacing  of  1  l.l  fs  and  a  pulse  length  of  0.21  fs. 

The  coherence  pab  (Fig.  1)  of  the  driven  molecular 
transition  is  central  to  this  work.  This  coherence  is 
established  by  driving  the  molecular  transition  slightly  off 
resonance  with  driving  lasers  of  sufficient  intensity  that 
the  product  of  their  Rabi  frequencies  exceeds  the  product 
of  the  detuning  from  the  molecular  electronic  states 
and  the  detuning  from  the  Raman  transition.  With  the 
linewidth  of  the  applied  laser  pulses  small  as  compared 
to  the  Raman  detuning,  the  magnitude  of  the  molecular 
coherence  approaches  0.5,  and  its  sign  is  determined  by 
the  sign  of  the  Raman  detuning. 

For  the  conditions  of  the  previous  paragraph,  the  gen¬ 
eration  and  phase-slip  lengths  are  of  the  same  order,  and 
Raman  sideband  generation  proceeds  collinearly  and  very 
differently  than  in  the  conventional  low  coherence  regime. 
In  essence,  the  molecular  motion  now  modulates  the  elec¬ 
tronic  refractive  index  much  in  the  same  way  that  a  low 
frequency  electric  field  modulates  the  refractive  index  of 
a  polar  crystal.  Incident  optical  frequencies  become  fre¬ 
quency  modulated  with  peak  (sinusoidal)  frequency  devia¬ 
tions  and  spectral  bandwidths  that  substantially  exceed  the 
width  of  the  visible  spectrum.  This  allows  the  extension 
of  compression  techniques,  such  as  the  use  of  group  ve¬ 
locity  dispersion  as  described  here,  to  a  new  regime  of 
short  pulse  compression. 

In  pertinent  prior  work,  (1)  the  authors  have  noted  how 
an  electromagnetically  induced  transparency  (EIT)-like 
excitation  may  be  used  to  generate  a  broad  comb  of  coher¬ 
ent  sidebands,  but  the  Bessel-function  nature  of  the  spec¬ 
trum  and  the  possibility  of  pulse  compression  were  not 


observed  [1].  (2)  Hakuta  and  colleagues  have  experi¬ 
mentally  demonstrated  the  generation  of  collinear  Raman 
sidebands  in  solid  molecular  hydrogen  [2].  (3)  Kawano 
et  cii  have  demonstrated  the  use  of  stimulated  rotational 
Raman  scattering  to  generate  a  broad  spectrum  and  dis¬ 
cussed  possibilities  for  mode  locking  this  spectrum  [3].  (4) 
There  has  been  considerable  work  in  the  area  of  on-axis 
Raman  sideband  generation  [4].  (5)  There  is  also  some 
relation  to  the  work  of  Kaplan  and  Shkolnikov  on  Raman 
solitons  [5],  (6)  Workers  in  the  field  of  high-order  har¬ 
monic  generation  have  noted  the  possibility  of  obtaining 
a  short  time  structure  [6].  (But,  in  high-order  harmonic 
generation,  less  than  one  part  in  10-4  of  the  incident  pulse 
energy  is  converted  into  the  phased  spectrum;  by  the  ap¬ 
proach  of  this  work,  near-complete  conversion  is  likely.) 
(7)  There  is  now  considerable  work  on  EIT  [7]  and  related 
subjects  [8].  (8)  It  has  been  shown  that,  by  using  a  large 
atomic  coherence,  nearly  complete  frequency  conversion 
will  occur  in  a  single  coherence  length  [9]. 

The  body  of  this  Letter  divides  into  four  parts:  (a)  a 
summary  of  the  formalism  describing  the  interaction  of 
the  multilevel  system  of  Fig.  1  with  a  comb  of  Raman 


FIG.  1.  Energy  level  schematic  for  establishing  coherence 
pab  in  a  molecular  system.  Laser  fields  are  applied  at  the 
frequencies  of  the  q  =  0  and  q  =  —  l  sidebands.  Here, 
the  detuning  from  Raman  resonance,  as  is  required  for 
establishing  the  antiphased  state,  is  negative. 
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sidebands  [1],  this  formalism  includes  all  dispersive 
effects:  (b)  the  FM  solution  for  the  ideal  case  of  negligible 
dispersion:  (c)  a  numerical  simulation  showing  concurrent 
spectral  generation  and  pulse  compression  in  Di\  (d)  a 
brief  discussion  of  measurement  possibilities. 

We  will  work  with  the  envelope  quantities  Eq(z,t)  and 
Pab(z.t)  for  electromagnetic  fields  and  density  matrix 
elements  such  that  the  full  time  and  space  quantities  are 
Eq(z,t)  =  Re{£9(z,f)  exp [j(coqt  -  k,z)]}  and  pfl*(z,t)  = 
Re{pfl*(z,r)  exp [j{comt  -  £mz)]},  with  coq  =  coo  + 
q(coi,  —  coa  —  A<u)  =  coo  +  qcom.  The  modulation 
frequency  com  or,  equivalently,  the  frequency  difference 
between  adjacent  sidebands,  is  equal  to  the  difference  of 
the  two  applied  frequencies  o»o  and  co~\.  The  two-photon 
detuning  A co  is  the  difference  of  com  and  the  Raman 
transition  frequency  and,  as  shown  in  Fig.  1,  is  negative; 


kq 


coq/c  and  km  =  co. 


,/c. 


We  allow  for  an  arbitrary  number  of  virtual  states  U) 
with  energies  R<at.  The  matrix  elements  from  states  |a) 
and  |fc)  to  these  states  are  p.ai  and  /i&(,  respectively.  When 
the  derivatives  of  the  probability  amplitudes  of  the  upper 
states  |i)  are  small  as  compared  to  the  detunings  from  these 
states,  the  problem  can  be  written  in  terms  of  an  effective, 
distance-dependent,  two-by-two  Hamiltonian  [1]: 


ffeff  = 


L  C  exp  (j^z) 


B  exp (-j-f-z) 
D  -  2A  co 


(1) 


,;C  =  B*;  and 
,(2.  All  rotating  and  nonrotating,  as  well  as 


where  A  =  Y.q  aq\Eq\2\  B  =  bqEqEq 
D  =  Z.qdq\E,‘ 
cross-transition,  terms  are  retained  within  this  Hamilton¬ 
ian.  We  also  assume  the  ideal  case  of  zero  linewidth  for 
the  |a)-|B)  transition.  The  constants  aq,  bq,  and  dq  deter¬ 
mine  the  dispersion  and  coupling  and  are 


aa  = 


ba  = 


1 

in2 


1 

in2 


i 

in2 


\^a 


(cOj  -  coa) 

I  Mml2 


CO  q 


( COj  -  coa)  +  t 

(COj  -  COa)  -  COq 
pai  jAbi 


'1  J 


+ 


(2) 


(cot  -  cob )  +  coq  J 

(l (Oi  -  COb)  -  (Oq 
\p-bi\2 


( (Oi  -  (Ob)  +  (Oq  J 

We  assume  that  all  of  the  molecular  population  is 
initially  in  the  ground  state  and  the  envelopes  of  the 
exciting  laser  pulses  vary  slowly  as  compared  to  the 
separation  of  the  eigenvalues  of  the  Hamiltonian  of  Eq.  ( 1 ). 
This  allows  adiabatic  preparation  of  all  of  the  molecules  in 
one  eigenstate.  By  choosing  the  sign  of  the  detuning  A  a>, 
either  the  phased  or  antiphased  eigenstate  is  selected  at  the 
input  of  the  medium  and  is,  as  shown  below,  maintained 
throughout  the  medium.  The  molecular  coherence  pab 


for  the  phased  state  has  the  same  sign  as  the  quantity  B 
in  Eq.  (1)  and,  for  the  antiphased  state,  has  the  opposite 
sign.  Both  eigenstates  lead  to  an  FM  signal  but,  as  shown 
below,  in  a  normally  dispersive  medium,  with  the  boundary 
condition  of  two  incident  frequencies,  only  the  antiphased 
state  allows  pulse  compression. 

Defining 

8  -  |fl|  e*pO>);  tan  0  “  2i„  o  +  A  ’  <3) 

the  eigenstate,  which  adiabatically  evolves  from  the 
ground  state  as  field  amplitudes  are  increased,  and  the 
coherence,  which  is  obtained  from  it,  are  as  follows: 

1+)  =  COS  y  exp^y 

+  sin  y  exp^-;'  yj  I b),  (4) 


Pab  =  y  sin  0  exp(y», 

where  positive  d  corresponds  to  the  phased  state  and 
negative  6  to  the  antiphased  state.  By  choosing  the 
intensities  of  the  applied  fields  and  the  detuning  Aw,  the 
quantity  sin  0  can  be  set  to  any  value  between  - 1  and  1 . 

The  propagation  equation  for  the  qth  sideband  in  local 
time  is 


dEq 

dz 


-jT)hcoqN(aqp 

aaEq 


"b  dqpbbEq  +  bqpabEq—  [ 


+  Cqp*abEq+ 1),  (5) 

where  N  is  the  number  of  molecules  per  volume  and 
V  =  Weo)l/2- 

In  order  to  understand  the  numerical  results  which 
follow,  we  first  exam  the  analytically  solvable  case  of 
negligible  dispersion  and  limited  modulation  bandwidth. 
We  take  all  of  the  sidebands  to  be  sufficiently  far  from 
the  resonances  that,  in  Eq,  (2),  aq  —  bq  =  cq  —  dq  = 
a o.  We  define  propagation  constants  fiq  =  r}fi(oqNao 
and  change  variables  by  Eq  =  Eq  Qxp(-j(3qz),  pab  = 
Pab  exp (-jfimz)-  With  these  approximations  and  defi¬ 
nitions,  Eq.  (5)  becomes 
dE 

=  —j  fiq(pabEq-\  +  pafyEq- fi).  (6) 


The  quantity  pab  =  \  sin  6  exp(y^),  where  9  and  <p 
are  given  by  Eq.  (3)  with  Eq  replaced  by  Eq.  As  the 
boundary  conditions  for  the  propagation  equation,  we 
assume  that  the  fields  £0(0)  and  £-i(0)  are  applied  at  the 
input  of  the  medium  and  all  other  Raman  sidebands  are 
generated  in  the  medium. 

Although  we  will  not  make  this  approximation  in 
the  numerical  work  which  follows,  here  we  assume 
a  restricted  modulation  bandwidth  such  that  all 


CO  Q 


are  equal  to  <o0  and  /3a  =  /3o.  By  using  the  Bessel- 
function  identities  l^fr-  =  [/n-iU)  ~  •Ai+iU)]  and 


Jn+q(x)Jn+P(x)  =  8pq,  it  may  be  verified  that,  with 
y  =  /3q  sin[0(O)],  the  solution  pab(z)  =  pab( 0)  and 
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Eq{z)  =  £o(0)  exp  j(^cp{ 0)  -  Jqiyz)  +*  £-i(0)  exp  ;^(0)  -  y)(<7  +  0  -Vi(yz)  (?) 

satisfies  both  density  matrix  and  propagation  equations  and  also  the  boundary  condition  at  z  —  0.  [<p(0)  is  obtained 
from  Eq.  (3).]  With  Jn{x)  exp (jncot)  =  exp[y*  sin(w/)],  the  time  domain  signal  is 

E(t)  =  Re{£'0(0)  exp  +  yz  cos((omt  +  <p(0))]  +  £-,(0)  exp  j[p-\t  +  yz  cos(o>mr  +  <p(0))]}.  (8) 


Equation  (8)  is  a  superposition  of  two  FM  signals  with 
center  frequencies  corresponding  to  the  frequencies  which 
are  applied  at  z  =  0  and  a  peak  phase  deviation  yz. 
It  is  important  that  the  frequency  modulation  leaves  the 
eigenvectors  of  the  Hamiltonian  unchanged. 

The  total  number  of  sidebands  generated  in  a  cell  of 
length  L  is  approximately  equal  to  2 y£.  The  assumption 
of  limited  modulation  bandwidth  (coq  =  gjo)  implies  that 
yz  «  (*)o/(Dm\  i.e.,  the  peak  frequency  excursion  must 
be  small  as  compared  to  the  center  frequency.  We 
have  checked  that,  when  this  assumption  is  satisfied,  the 
solution  of  Eq.  (7)  coincides  with  the  exact  numerical 
solution  of  Eq.  (6). 

Equations  (7)  and  (8)  hold  for  both  the  antiphased 
(sin  0  <  0)  and  the  phased  case  (sin  0  >  0).  But,  in  the 
antiphased  case,  the  higher  frequencies  (blue)  are  tempo¬ 
rally  advanced  with  regard  to  the  lower  frequencies  (red), 
while,  in  the  phased  case,  red  is  advanced  as  compared 
to  blue.  When  £0(0)  -  E-i(0),  only  the  signal  gener¬ 
ated  by  the  antiphased  state  allows  pulse  compression  in  a 
normally  dispersive  medium.  Figure  2(a)  shows  the  time 
domain  signal  produced  by  the  antiphased  state;  i.e.,  an 
FM  signal  with  an  envelope  at  the  beat  frequency  (om . 

Pulse  compression  is  obtained  by  setting  the  group  ve¬ 
locity  delay  which  corresponds  to  the  frequency  deviation 


time  (i/<om) 

FIG.  2.  Electric  field  amplitude  vs  time  (a)  as  produced  by 
an  antiphased  state  without  pulse  compression.  The  parame¬ 
ters  are  £q(0)  =  £-i(0)  =  1,  a>o/a)m  =  18,  and  yL  =  7. 
(b)  Following  pulse  compression  in  a  medium  with  a 
frequency-invariant  group  velocity  dispersion  and  a  length 
determined  by  Eq.  (9). 


cod  —  wmyL  equal  to  Tr/2a)m.  Here,  but  not  in  the  nu¬ 
merical  work  which  follows,  we  assume  a  frequency  in¬ 
variant  group  velocity  dispersion  d2k/da)2.  The  medium 
length  /  which  is  necessary  for  pulse  compression  is  then 


Figure  2(b)  shows  the  signal  which  is  obtained  when  the 
signal  of  Fig.  2(a)  is  transmitted  through  a  medium  of  this 
length.  Pulse  compression  with  a  ratio  of  pulse  width  to 
pulse  separation  of  about  1/2 yL  is  observed. 

We  turn  next  to  the  numerical  solution  for  D2.  We  do 
not  make  any  of  the  dispersive  approximations  described 
above.  We  find  that  both  generation  and  pulse  com¬ 
pression  may  be  accomplished  concurrently  in  the  same 
medium.  We  assume  Eq  and  £-i  to  be  applied  at  z  =  0 
and  retain  sidebands  between  q  —  -9  and  (7  =  18.  The 
solution  of  Eq.  (5)  is  obtained  by  forward  stepping  from 
z  =  0  with  the  density  matrix  elements  recalculated  at 
each  step. 

We  consider  the  fundamental  vibrational  transition  in 
D2  with  <om  =  2994  cm"1  and  take  all  molecules  in  the 
7  =  0  ground  state.  The  constants  aq,bq,  and  dq  include 
the  contributions  of  the  0  -  52nd  vibrational  transitions 
of  the  Lyman  band  and  the  0  -  18th  transitions  of  the 
Werner  band.  Oscillator  strengths  and  Franck-Condon 
factors  are  obtained  from  Allison  and  Dalgamo  and 
energies  are  obtained  from  Herzberg  [10].  The  £>2  density 
and  the  total  cell  length  are  chosen  such  that  the  zeroth 
sideband,  if  propagating  alone,  would  accumulate  a  phase 
shift  fio L  =  147r  rad.  At  a  pressure  of  1  atm,  this 
corresponds  to  L  =  3.6  cm. 

We  envision  using  single-mode  pulsed  lasers  with  pulse 
lengths  of  —5  ns  and  pulse  energies  of  -50  mJ.  We 
take  the  applied  laser  frequencies  to  be  <o  =  28  183  cm-1 
(frequency-tripled  Nd:YAG)  and  o)~\  =  25  189  cm  1 
(frequency-doubled  Ti:sapphire).  We  choose  the  Raman 
detuning  =  - 1  GHz  to  be  much  larger  than  both  the 
Raman  linewidth  (~0.3  GHz  for  D2  at  STP)  and  the  laser 
linewidth  (—0.1  GHz  for  a  single-mode,  5-ns-pulse  laser). 
We  set  the  power  densities  of  £0  and  £-i  at  z  =  0  to 
10u  W/cm2.  Evaluating  Eqs.  (l)-(4),  sin  0(0)  »  -0.9. 
The  coherence  \pab\  ~  0.45  and  the  cell  length  are  set 
such  that  about  twenty  sidebands  are  generated,  and 
the  concurrent  group  velocity  delay  is  correct  for  pulse 


compression. 

Figure  3  shows  spectral  and  temporal  evolution  as  a 
function  of  distance  for  the  conditions  of  the  previous 
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frequency  (1,000  cm"')  time  (fs) 


FIG.  3.  Frequency  spectrum  (left)  and  instantaneous  power 
density  vs  time  (right)  in  D2  at  1  atm  pressure  with  applied 
power  densities  at  z  =  0  such  that  \pab\  “  0-45.  Parts  (a)-(c) 
show  the  spectrum  and  temporal  waveform  for  (a)  z  =  0  cm, 
(b)  z  =  1.8  cm,  and  (c)  :  =  3.6  cm.  Sideband  powers  are 
normalized  to  those  of  the  incident  fields.  The  dashed  lines  on 
the  time  domain  plots  show  the  envelopes  of  possible  single¬ 
cycle  pulses. 


paragraph.  Sideband  frequencies  are  predicted  to  range 
from  about  1.4  p. m  in  the  infrared  to  150  nm  in  the 
vacuum  ultraviolet.  At  the  appropriate  cell  length  these 
frequencies  synthesize  a  train  of  pulses  with  a  pulse  sepa¬ 
ration  of  11.1  fs  and  a  pulse  length  of  0.21  fs.  We  find 
sharp  pulse  compression  at  experimental  conditions  which 
are  readily  realizable.  Because  the  ratio  of  the  applied 
optical  frequencies  to  the  Raman  transition  frequency  is, 
in  general,  noninteger,  the  peak  of  the  optical  cycle  will 
drift  from  pulse  to  pulse  within  a  narrow  envelope.  This 
envelope  (width  =  0.72  fs)  is  shown  by  the  dashed  line 
in  each  of  the  time  domain  plots  of  Fig.  3.  For  the 
numerical  analysis  we  have  treated  the  driving  fields  as 
monochromatic.  For  the  pulsed  excitation,  as  described 
above,  we  expect  these  results  to  apply  for  a  time  about 
equal  to  the  dephasing  time  of  the  Raman  transition,  and 
therefore  expect  nearly  all  of  the  incident  power  to  be  con¬ 
verted  into  the  periodic  train  of  subfemtosecond  pulses. 

The  predictions  of  this  numerical  analysis  should  be 
verifiable  by  two  methods.  The  relative  intensities  of  the 
sidebands  could  be  measured;  these  differ  in  the  phased 
and  antiphased  states,  and  are  correlated  with  the  phase 
of  each  sideband.  In  the  time  domain,  the  peak  power 
density  grows  on  a  scale  of  3  cm  and  varies  in  a  predictable 
manner  thereafter.  In  the  region  of  the  highest  peak  power 
(~5  mm  long),  one  might  expect  multiphoton  processes  in 
either  D2  itself  or  in  a  dopant  to  light  up  as  an  indication 


of  the  presence  of  ultrashort  pulses.  The  physics  of  such 
processes,  as  well  as  other  nonlinear  behavior  in  this 
regime  is,  in  its  own  right,  of  considerable  interest. 

The  principle  approximations  of  the  numerical  results 
given  here  are  the  neglect  of  the  dephasing  of  the 
Raman  transition,  the  assumption  that  the  driving  fields 
are  both  quasimonochromatic  and  of  an  infinite  extent 
in  the  transverse  direction,  and  that  all  molecules  are 
initially  in  the  J  =  0  ground  state.  Although  the  effects 
of  self-focusing,  temporal  chirping,  saturation  of  the 
Raman  transition,  and  thermal  population  have  not  been 
discussed,  we  do  not  see  contradictions. 

In  summary,  we  have  described  how  the  ideas  of  fre¬ 
quency  modulation  and  pulse  compression  may  be  ex¬ 
tended  to  the  atomic  time  scale.  The  technique  will 
allow  the  generation  of  pulses  of  high  intensity  radiation 
with  a  pulse  length  which  is  less  than  the  revolution  time 
of  an  electron  in  the  first  Bohr  orbit  of  atomic  hydrogen. 
The  use  of  rotational  Raman  transitions,  instead  of  vibra¬ 
tional  transitions,  is  promising.  For  example,  the  v'  =  0, 
j1  =  2  — 1 ►  v"  —  0,  J"  =  0  rotational  transition  in  H2  will 
produce  subfemtosecond  pulses  with  a  pulse  spacing  of 
94  fs  [11]. 

This  work  was  supported  by  the  U.S.  Army  Research 
Office,  the  U.S.  Air  Force  Office  of  Scientific  Research, 
and  the  U.S.  Office  of  Naval  Research. 
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Techniques  that  use  quantum  interference  effects  are  being 
actively  investigated  to  manipulate  the  optical  properties  of 
quantum  systems1.  One  such  example  is  electromagnetically 
induced  transparency,  a  quantum  effect  that  permits  the  propaga¬ 
tion  of  light  pulses  through  an  otherwise  opaque  medium2"5.  Here 
we  report  an  experimental  demonstration  of  electromagnetically 
induced  transparency  in  an  ultracold  gas  of  sodium  atoms,  in 
which  the  optical  pulses  propagate  at  twenty  million  times  slower 
than  the  speed  of  light  in  a  vacuum.  The  gas  is  cooled  to 
nanokelvin  temperatures  by  laser  and  evaporative  cooling6"10. 
The  quantum  interference  controlling  the  optical  properties  of 
the  medium  is  set  up  by  a  ‘coupling’  laser  beam  propagating  at  a 
right  angle  to  the  pulsed  ‘probe’  beam.  At  nanokelvin  tempera¬ 
tures,  the  variation  of  refractive  index  with  probe  frequency  can 
be  made  very  steep.  In  conjunction  with  the  high  atomic  density, 
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this  results  in  the  exceptionally  low  light  speeds  observed.  By 
cooling  the  cloud  below  the  transition  temperature  for  Bose- 
Emstein  condensation11" 15  (causing  a  macroscopic  population  of 
alkali  atoms  in  the  quantum  ground  state  of  the  confining 
potential),  we  observe  even  lower  pulse  propagation  velocities 
(17ms1)  owing  to  the  increased  atom  density.  We  report  an 
inferred  nonlinear  refractive  index  of  0.18  cm2  W"1  and  find  that 
the  system  shows  exceptionally  large  optical  non  linearities,  which 
are  of  potential  fundamental  and  technological  interest  for  quantum 
optics. 

The  experiment  is  performed  with  a  gas  of  sodium  atoms  cooled 
to  nanokelvin  temperatures.  Our  atom  cooling  set-up  is  described 
in  some  detail  in  ref.  14.  Atoms  emitted  from  a  candlestick’  atomic 
beam  source15  are  decelerated  in  a  Zeeman  slower  and  loaded  into  a 
magneto-optical  trap.  In  a  few  seconds  we  collect  a  cloud  of  1010 
atoms  at  a  temperature  of  1  mK  and  a  density  of  6  X  10n  cm-3.  The 
atoms  are  then  polarization  gradient  cooled  for  a  few  milliseconds 
to  50  p.K  and  optically  pumped  into  the  F  =  1  ground  state  with  an 
equal  population  of  the  three  magnetic  sublevels.  We  then  turn  all 
laser  beams  off  and  coniine  the  atoms  magnetically  in  the  '4  Dee’ 
trap14.  Only  atoms  in  the  MF  =  -  1  state,  with  magnetic  dipole 
moments  directed  opposite  to  the  magnetic  field  direction  (picked 
as  the  quantization  axis),  are  trapped  in  the  asymmetric  harmonic 
trapping  potential.  This  magnetic  filtering  results  in  a  sample  of 
atoms  that  are  all  in  a  single  atomic  state  (state  |  1)  in  Fig.  lb)  which 
allows  adiabatic  optical  preparation  of  the  atoms,  as  described 
below,  and  minimal  heating  of  the  cloud. 

Next  we  evaporativelv  cool  the  atoms  for  38  s  to  the  transition 
temperature  for  Bose-Einstein  condensation,  Tc.  The  magnetic 
fields  are  then  adjusted  to  adiabatically  soften  the  trap.  The  resulting 
trapping  potential  has  a  frequency  of  /.  =  21  Hz  along  the  sym¬ 
metry  ( z )  axis  of  the  4  Dee  trap,  and  transverse  frequencies 
f  x  =  /,  =  69  Hz.  The  bias  field,  parallel  to  the  z  axis,  is  11G. 
When  we  cool  well  below  Tc,  we  are  left  with  1-2  million  atoms 
in  the  condensate.  For  these  parameters  the  transition  occurs  at  a 
temperature  of  Ic  =  435  nK  and  a  peak  density  in  the  cloud  of 
5  X  1012 cm"3. 

We  now  apply  a  linearly  polarized  laser  beam,  the  coupling  beam, 
tuned  to  the  transition  between  the  unpopulated  hyperfine  states  |  2) 
and  |  3)  (Fig.  lb).  This  beam  couples  states  |  2)  and  |  3)  and  creates  a 
quantum  interference  for  a  weaker  probe  laser  beam  (left  circularly 
polarized)  which  is  tuned  to  the  1 1)  — ►  |3)  transition.  A  stable 
eigenstate  (the  ’dark  state’)  of  the  atom  in  the  presence  of  coupling 
and  probe  lasers  is  a  coherent  superposition  of  the  two  hyperfine 
ground  states  |  l)  and  |  2).  The  ratio  of  the  probability  amplitudes  is 
such  that  the  contributions  to  the  atomic  dipole  moment  induced 
by  the  two  lasers  exactly  cancel.  The  quantum  interference  occurs  in 
a  narrow  interval  of  probe  frequencies,  with  a  width  determined  by 
the  coupling  laser  power. 

Figure  2a  shows  the  calculated  transmission  of  the  probe  beam  as 
a  function  of  its  detuning  from  resonance  for  parameters  which  are 
typical  of  this  work.  In  the  absence  of  dephasing  of  the  |l)— *■  |2) 
transition,  the  quantum  interference  would  be  perfect,  and  at  line 
centre,  the  transmission  would  be  unity.  Figure  2b  shows  the 
refractive  index  for  the  probe  beam  as  a  function  of  detuning. 
Due  to  the  very  small  Doppler  broadening  of  the  1 1)— » ►  1 2)  transition 
in  our  nanokelvin  samples,  application  of  very  low  coupling 
intensity  leads  to  a  transparency  peak  with  a  width  much  smaller 
than  the  natural  line  width  of  the  j  1)  — ►  |3)  transition.  Correspond¬ 
ingly,  the  dispersion  curve  is  much  steeper  than  can  be  obtained  by 
any  other  technique,  and  this  results  in  the  unprecedented  low 
group  velocities  reported  here.  The  group  velocity  for  a  propagat¬ 
ing  electromagnetic  pulse  is16-19: 
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Figure  1  Experimental  set-up.  A  'coupling'  laser  beam  propagates  along  the  x 
axis  with  its  linear  polarization  along  the  11-G  bias  field  in  thez  direction.  The 
probe'  laser  pulse  propagates  along  thez  axis  and  is  ieft-circuiady  polarized. 
With  a  flipper  mirror  in  front  of  the  camera  CCD  1,  we  direct  this  probe  beam  either 
to  the  camera  or  to  the  photomultiplier  (PMT).  For  pulse  delay  measurements,  we 
oiace  a  pinhole  in  an  external  image  plane  of  the  imaging  optics  and  select  a  small 
area,  iB^m  in  diameter,  of  the  probe  beam  centred  on  the  atom  clouds  (as 
indicated  by  the  dashed  circle  in  inset  (i)).  The  pulse  delays  are  measured  with  the 
PMT.  The  imaging  beam  propagating  along  the  y  axis  is  used  to  image  atom 
ciouos  onto  camera  CCD  2  to  find  the  length  of  the  clouds  along  the  pulse 
propagation  direction  (z  axis)  for  determination  of  light  speeds,  inset  (ii)  shows 
atoms  cooied  to  450  nK  which  is  i5nK  above  Tr,  (Note  that  this  imaging  beam  is 
never  applied  at  the  same  time  as  the  probe  pulse  and  coupling  laser).  The 
position  of  a  cloud  and  its  diameter  in  the  two  transverse  directions, x  andy,  are 
found  with  CCD  1.  Inset  (i)  shows  an  image  of  a  condensate. 
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:  Hfcre  «(o)p)  is  the  refractive  index  at  probe  frequency  o>p  (rads-1), 
|.QC|2  is  the  square  of  the  Rabi  frequency  for  the  coupling  laser  and 
varies  linearly  with  intensity,  |x13  is  the  electric  dipole  matrix 
|  element  between  states  |  1)  and  |  3),  N  is  the  atomic  density,  and 
|  e0  is  the  permittivity  of  free  space.  At  line  centre,  the  refractive  index 
;  is  unity,  and  the  second  term  in  the  denominator  of  equation  (1) 
t  dominates  the  first.  An  important  characteristic  of  the  refractive 
;  index  profile  is  that  on  resonance  the  dispersion  of  the  group 
|  velocity  is  zero  (see  ref.  16),  that  is,  d2n/dco2  =  0,  and  to  lowest 
;  order,  the  pulse  maintains  its  shape  as  it  propagates.  The  established 
|  quantum  interference  allows  pulse  transmission  through  our  atom 
clouds  which  would  otherwise  have  transmission  coefficients  of  e_n0 
(below  rc)>  and  creates  a  steep  dispersive  profile  and  very  low  group 
velocity  for  light  pulses  propagating  through  the  clouds. 

We  note  that  the  centres  of  the  curves  in  Fig.  2  are  shifted  by 
0.6  MHz  from  probe  resonance.  This  is  due  to  a  coupling  of  state  |  2) 
to  state  |  4)  through  the  coupling  laser  field,  which  results  in  an  a.c. 
Stark  shift  of  level  |  2)  and  a  corresponding  line  shift  of  the  2  — ►  3 
transition.  As  the  transparency  peak  and  unity  refractive  index  are 
obtained  at  two-photon  resonance,  this  leads  to  a  refractive  index  at 
j  the  1  — +  3  resonance  frequency  which  is  different  from  unity.  The 
difference  is  proportional  to  the  a.c.  Stark  shift  and  hence  to  the 
coupling  laser  intensity,  which  is  important  for  predicting  the 
nonlinear  refractive  index  as  described  below. 

A  diagram  of  the  experiment  is  shown  in  Fig.  la.  The  2.5-mm- 
diameter  coupling  beam  propagates  along  the  x  axis  with  its  linear 
;  polarization  parallel  to  the  B  field.  The  0.5-mm-diameter,  a~ 
polarized  probe  beam  propagates  along  the  z  axis.  The  size  and 
position  of  the  atom  cloud  in  the  transverse  directions,  x  and  y,  are 
obtained  by  imaging  the  transmission  profile  of  the  probe  beam 
after  the  cloud  onto  a  charge-coupled-device  (CCD)  camera.  An 
image  of  a  condensate  is  shown  as  inset  (i).  A  55  mW  cm'2  coupling 


laser  beam  was  present  during  the  10-p.s  exposure  of  the  atoms  to  a  \ 
5  mW  cm”2  probe  beam  tuned  close  to  resonance.  The//7  imaging  j 
optics  are  diffraction -limited  to  a  resolution  of  7  jxm.  ! 

During  the  pulse  delay  experiments,  a  pinhole  (placed  in  an 
external  image  plane  of  the  lens  system)  is  used  to  select  only  the 
part  of  the  probe  light  that  has  passed  through  the  central  15  of 
the  atom  cloud  where  the  column  density  is  the  greatest.  The  outline 
of  the  pinhole  is  indicated  with  the  dashed  circle  in  inset  (i). 

Both  coupling  and  probe  beams  are  derived  from  the  same  dye 
laser.  The  frequency  of  the  coupling  beam  is  set  by  an  acousto-optic 
modulator  (AOM)  to  the  |2)  — ►  |3)  resonance.  Here  we  take  into 
account  both  Zeeman  shifts  and  the  a.c.  Stark  shift  described  above. 

The  corresponding  probe  resonance  is  found  by  measuring  the 
transmission  of  the  probe  beam  as  a  function  of  its  frequency.  We 
apply  a  fast  frequency  sweep,  across  32  MHz  in  50  \ xs,  and  determine 
resonance  from  the  transmission  peak.  The  sweep  is  controlled  by  a 
separate  AOM.  The  frequency  is  then  fixed  at  resonance,  and  the 
temporal  shape  of  the  probe  pulse  is  generated  by  controlling  the  r.f.  ! 
drive  power  to  the  AOM.  The  resulting  pulse  is  approximately  j 
gaussian  with  a  full-width  at  half-maximum  of  2.5  \xs.  The  peak  j 
power  is  1  mW  cm'2  corresponding  to  a  Rabi  frequency  of 
Op  =  0.20  A,  where  the  Einstein  A  coefficient  is  6.3  X  107  rads-1,  | 
To  avoid  distortion  of  the  pulse,  it  is  made  of  sufficient  duration  that 
its  Fourier  components  are  contained  within  the  transparency  peak. 

Probe  pulses  are  launched  along  the  z  axis  4  |xs  after  the  coupling  j 
beam  is  turned  on  (the  coupling  field  is  left  on  for  100  |is).  Due  to  I 
the  magnetic  filtering  discussed  above,  all  atoms  are  initially  in  state  | 
|  1)  which  is  a  dark  state  in  the  presence  of  the  coupling  laser  only. 
When  the  pulse  arrives,  the  atoms  adiabatically  evolve  so  that 
the  probability  amplitude  of  state  |  2)  is  equal  to  the  ratio 
f2p/(f2p  +  Q\)m,  where  Qp  is  the  probe  Rabi  frequency.  To  establish 
the  coherent  superposition  state,  energy  is  transferred  from  the 


|  Figure  2  Effect  of  probe  detuning,  a,  Transmission  profile.  Calculated  probe 
;  transmission  as  a  function  of  detuning  from  the  |1>—  |3)  resonance  for  an  atom 
i  cloud  cooled  to  450 nK,  with  a  peak  density  of  3.3  x  1012  cm _3  and  a  length  of 
229  n-m  (corresponding  to  the  cloud  in  inset  (ii)  of  Fig.  ta).  The  coupling  laser  is 
|  resonant  with  the  |2>— 13)  transition  and  has  a  power  density  of  52  mWcm-2.  b, 
j  Refractive  index  profile.  The  calculated  refractive  index  is  shown  as  a  function  of 
1  probe  detuning  for  the  same  parameters  as  in  a.  The  steepness  of  the  slope  at 
.  resonance  is  inversely  proportional  to  the  group  velocity  of  transmitted  light 
j  pulses  and  is  controlled  by  the  coupling  laser  intensity.  Note  that  as  a  result  of  the 
:  a.c.  Stark  shift  of  the  |2>  —  |3)  transition,  caused  by  a  coupling  of  states  1 2)  and  1 4> 
;  through  the  coupling  laser  field,  the  centre  of  the  transmission  and  refractive 
:  index  profiles  is  shifted  by  0.6  MHz.  The  shift  of  the  refractive  index  profile  results 
i  in  the  nonlinear  refractive  index  described  in  the  text. 


i 

I 


Figure  3  Pulse  delay  measurement.  The  front  pulse  (open  circles)  is  a  reference 
pulse  with  no  atoms  in  the  system.  The  other  pulse  (filled  circles)  is  delayed  by 
7.05  m.s  in  a  229-fim-long  atom  cloud  (see  inset  (ii)  in  Fig.  la).  The  corresponding 
light  speed  is  32.5 ms"'.  The  curves  represent  gaussian  fits  to  the  measured 
pulses. 
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fronj  °f  fhe  probe  pulse  to  the  atoms  and  the  coupling  laser  field.  At 
the  end  of  the  pulse,  the  atoms  adiabatically  return  to  the  original 
state  |  l)  and  the  energy  returns  to  the  back  of  the  probe  pulse  with 
no  net  energy  and  momentum  transfer  to  the  atomic  cloud.  Because 
the  refractive  index  is  unity,  the  electric  field  is  unchanged  as  the 
probe  pulse  enters  the  medium.  As  the  group  velocity  is  decreased, 
the  total  energy  density  must  increase  so  as  to  keep  constant  the 
power  per  area.  This  increase  is  represented  by  the  energy  stored  in 
,  the  atoms  and  the  coupling  laser  field  during  pulse  propagation 
'  through  the  cloud. 

The  pulses  are  recorded  with  a  photomultiplier  (3-ns  response 
time)  after  they  penetrate  the  atom  clouds.  The  output  from  the 
photomultiplier  is  amplified  by  a  150-MHz-bandwidth  amplifier 
and  the  waveforms  are  recorded  on  a  digital  scope.  With  a  ‘flipper 
mirror  in  front  of  the  camera  we  control  whether  the  probe  beam  is 
directed  to  the  camera  or  to  the  photomultiplier. 

The  result  of  a  pulse  delay  measurement  is  shown  in  Fig.  3.  The 
front  pulse  is  a  reference  pulse  obtained  with  no  atoms  present.  The 
pulse  delayed  by  7.05  fis  was  slowed  down  in  an  atom  cloud  with  a 
length  of  229  |xm  (see  Fig.  la,  inset  ( ii) ) .  The  resulting  light  speed  is 
32.5  ms”1.  We  used  a  coupling  laser  intensity  of  12mWcm“2 
corresponding  to  a  Rabi  frequency  of  Oc  =  0.56  A.  The  cloud  was 
cooled  to  450  nK  (which  is  15  nK  above  Tc),  the  peak  density  was 
3.3  X  10l2cm”\  and  the  total  number  of  atoms  was  3.8  X  106. 
From  these  numbers  we  calculate  that  the  pulse  transmission 
coefficient  would  be  e”63  in  the  absence  of  the  coupling  laser.  The 
probe  pulse  was  indeed  observed  to  be  totally  absorbed  by  the  atoms 
when  the  coupling  beam  was  left  off.  Inhomogeneous  broadening 
due  to  spatially  varying  Zeeman  shifts  is  negligible  (^20 kHz)  for 
the  low  temperatures  and  correspondingly  small  cloud  sizes  used 
here. 

The  size  of  the  atom  cloud  in  the  z  direction  is  obtained  with 
another  CCD  camera.  For  this  purpose,  we  use  a  separate  1  mW  cm”2 
laser  beam  propagating  along  the  vertical  y  axis  and  tuned  20  MHz 
below  the  F  =  2  — » ►  3  transition.  The  atoms  are  pumped  to  the 
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Figure  4  Light  speed  versus  atom  cloud  temperature.  The  speed  decreases  with 
temperature  due  to  the  atom  density  increase.  The  open  circles  are  for  a  coupling 
power  of  52  mW  cm"2  and  the  filled  circles  are  for  a  coupling  power  of  12  mW  cm"2. 
The  temperature  7L  marks  the  transition  temperature  for  Bose-Einstein 
condensation.  The  decrease  in  group  velocity  below  Tc  is  due  to  a  density 
increase  of  the  atom  cloud  when  the  condensate  is  formed.  From  imaging 
measurements  we  obtain  a  maximum  atom  density  of  8  x  I0,3cm-3  at  a  tem¬ 
perature  of  200  nK.  Here,  the  dense  condensate  component  constitutes  60%  of  ail 
atoms,  and  the  total  atom  density  is  16  times  larger  than  the  density  of  a  non- 
condensed  cloud  at  rc.  The  light  speed  measurement  at  50  nK  is  for  a  cloud  with  a 
condensate  fraction  5=90%.  The  finite  dephasing  rate  due  to  state  |  4}  does  not 
allow  pulse  penetration  of  the  most  dense  clouds.  This  problem  could  be  over¬ 
come  by  tuning  the  laser  to  the  D,  line  as  described  in  the  text. 


F  —  2  ground  state  for  10  fxs  before  the  imaging  which  is  performed  j 
with  an  exposure  time  of  10  pis.  We  image  the  transmission  profile  | 
of  the  laser  beam  after  the  atom  cloud  with  diffraction-limited //5  j 
optics.  An  example  is  shown  in  Fig.  la,  inset  (ii),  where  the  j 
asymmetry  of  the  trap  is  clear  from  the  cloud’s  elliptical  profile. 
We  note  that  the  imaging  laser  is  never  applied  at  the  same  time  as 
the  coupling  laser  and  probe  pulse,  and  for  each  recorded  pulse  or  j 
CCD  picture  a  new  cloud  is  loaded. 

We  measured  a  series  of  pulse  delays  and  corresponding  cloud 
sizes  for  atoms  cooled  to  temperatures  between  2.5  jxK  and  50  nK.  \ 
From  these  pairs  of  numbers  we  obtain  the  corresponding  propaga-  ! 
tion  velocities  (Fig.  4).  The  open  circles  are  for  a  coupling  power  of  j 
52  mW  cm  2  (f2c  =  1 .2  A).  The  light  speed  is  inversely  proportional  j 
to  the  atom  density  (equation  (1))  which  increases  with  lower  j 
temperatures,  with  an  additional  density  increase  when  a  conden-  ! 
sate  is  formed.  The  filled  circles  are  for  a  coupling  power  of  j 
12  mW  cm-2.  The  lower  coupling  power  is  seen  to  cause  a  decrease  ! 
of  group  velocities  in  agreement  with  equation  ( 1 ) .  We  obtain  a  light 
speed  of  17  ms”1  for  pulse  propagation  in  an  atom  cloud  initially 
prepared  as  an  almost  pure  Bose-Einstein  condensate  (condensate  i 
fraction  is  2*90%).  Whether  the  cloud  remains  a  condensate  during  , 
and  after  pulse  propagation  is  an  issue  that  is  beyond  the  scope  of 
this  Letter. 

Transitions  from  state  |  2)  to  state  |  4),  induced  by  the  coupling 
laser  (detuned  by  60  MHz  from  this  transition),  result  in  a  finite 
decay  rate  of  the  established  coherence  between  states  |  i)  and  |  2) 
and  limit  pulse  transmission.  The  dephasing  rate  is  proportional  to 
the  power  density  of  the  coupling  laser  and  we  expect,  and  find,  that 
probe  pulses  have  a  peak  transmission  that  is  independent  of 
coupling  intensity  and  a  velocity  which  reduces  linearly  with  this 
intensity.  The  dephasing  time  is  determined  from  the  slope  of  a 
semi-log  plot  of  transmission  versus  pulse  delay19.  At  a  coupling 
power  of  12  mW  cm”2,  we  measured  a  dephasing  time  of  9  juls  for  ! 
atom  clouds  just  above  Tc, 

Giant  Kerr  nonlinearities  are  of  interest  for  areas  of  quantum  j 
optics  such  as  optical  squeezing,  quantum  nondemolition,  and  ; 
studies  of  nonlocality.  It  was  recently  proposed  that  they  may  be  1 
obtained  using  electromagnetically  induced  transparency20.  Here  we  j 
report  the  first  (to  our  knowledge)  measurement  of  such  a  non-  j 
linearity.  The  refractive  index  for  zero  probe  detuning  is  given  by  i 
n  —  1  +  (n2Jc)  where  Ic  is  the  coupling  laser  intensity,  and  n2  the  i 
cross  phase  nonlinear  refractive  index.  As  seen  from  Fig.  2b,  the  : 
nonlinear  term  (n2Ic)  equals  the  product  of  the  slope  of  the 
refractive  index  at  probe  resonance  and  the  a.c.  Stark  shift  of  the  ! 
|2) — ►  |3>  transition  caused  by  the  coupling  laser.  We  can  then 
express  n2  by  the  formula  (see  equation  (1)); 

A(j)s  d  n  1  Auc  X  ! 

»2  =  -TT-  -  r— r1-  (2)  I 

lc  dcop  2x  7C  vg 

where  Aojs  is  the  a.c.  Stark  shift,  proportional  to  Ia  and  X  the 
wavelength  of  the  probe  transition.  We  measured  an  a.c.  Stark  shift  j 
of  1.3  X  106  rad  s  1  for  a  coupling  laser  intensity  of  40  mW  cm”2. 
For  a  measured  group  velocity  of  17ms“l  (Fig.  4),  we  obtain  a 
nonlinear  refractive  index  of  0.18  cm2  W~l.  This  nonlinear  index  is 
— 106  times  greater  than  that  measured  in  cold  Cs  atoms21.  j 

With  a  system  that  avoids  the  |  l)-|  2)  dephasing  rate  described  ! 
above  (which  can  be  obtained  by  tuning  to  the  D\  line  in  sodium), 
the  method  used  here  could  be  developed  to  yield  the  collision- 
induced  dephasing  rate  of  the  double  condensate  which  is  generated  | 
in  the  process  of  establishing  electromagnetically  induced  transpar¬ 
ency  (see  also  refs  22,  23).  In  that  case,  the  square  of  the  probability 
amplitude  for  state  |  3)  could  be  kept  below  10”3  during  pulse 
propagation,  with  no  heating  of  the  condensate  as  a  result.  With  I 
improved  frequency  stability  of  our  set-up  and  lower  coupling  : 
intensities,  even  lower  light  speeds  would  be  possible,  perhaps  of 
the  order  of  centimetres  per  second,  comparable  to  the  speed  of 
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jSound  in  a  Bose-Einstein  condensate.  Under  these  conditions  we 
expect  phonon  excitation  during  light  pulse  propagation  through 
the  condensate.  By  deliberately  tuning  another  laser  beam  to  the 
|2)— *•  1 4)  transition,  it  should  be  possible  to  demonstrate  optical 
switching  at  the  single  photon  level24.  Finally,  we  note  that  during 
propagation  of  the  atom  clouds,  light  pulses  are  compressed  in  the  z 
direction  by  a  ratio  of  clvg.  For  our  experimental  parameters,  that 
results  in  pulses  with  a  spatial  extent  of  only  43  p.m.  □ 
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We  demonstrate  efficient,  pulsed,  gas-phase,  nonlinear  frequency  conversion  in  a  quadrupiy  resonant, 
double-A  system  and,  simultaneously,  verify  theoretical  predictions  of  Rabi-frequency  matching  unique 
to  absorbing  nonlinear  media.  This  system  is  used  to  up-convert  ultraviolet  light  at  233  nm  to  the 
vacuum  ultraviolet  at  186  nm  in  atomic  Pb  vapor  with  small-signal  conversion  efficiencies  exceeding 
30%  and  with  modest  atomic  density-length  (NL)  products  (scale  1014  cm-2)  and  optical  power  densities 
(10-100  kW/cm2). 

PACS  numbers:  42.50.Gy,  42.50.Hz.  42.65.Ky,  42.79.Nv 


It  has  been  shown  that  the  techniques  of  electromagnet- 
ically  induced  transparency  (EIT)  may  greatly  improve  the 
efficiency  of  gas-phase  nonlinear  optical  processes  [1,2]. 
EIT  facilitates  the  preparation  of  near-maximal  coherence 
of  a  dipole-forbidden  transition  without  absorption  or  self¬ 
phase  modulation  of  the  applied  fields;  the  phase-coherent 
atomic  ensemble  thus  created  acts  as  a  strong  local  oscil¬ 
lator  which  may  mix  with  the  spectral  components  of  the 
applied  fields  to  generate  a  corresponding  spectrum  of  sum 
or  difference  frequencies.  Previous  work  in  both  the  gas 
[1,2]  and  solid  [3]  phases  has  focused  on  the  reactive  case 
(i.e.,  large  detuning  of  the  generated  frequency  from  any 
resonances  and  near-lossless  propagation)  and  has  demon¬ 
strated  efficient  parametric  frequency  conversion  to  the  ul¬ 
traviolet  and  vacuum  ultraviolet  (VUV). 

In  this  Letter,  we  explore  gas-phase,  nonlinear  frequency 
up-conversion  in  a  quadrupiy  resonant  atomic  system.  This 
atom-field  combination  is  known  as  a  resonant  double-A 
configuration  and  has  been  studied  in  the  context  of  EIT 
[4],  lasers  without  inversion  [5],  and  nonlinear  optics  [6]. 
The  short-wavelength  up-converted  radiation  is  generated 
directly  on  resonance  under  conditions  where,  in  the  ab¬ 
sence  of  the  other  fields,  it  would  be  strongly  absorbed. 
We  demonstrate,  however,  that  through  the  up-conversion 
process,  the  atomic  system  renders  itself  transparent  to  all 
four  fields.  Transparency  is  achieved  once  the  complex 
ratio  of  Rabi  frequencies  in  each  A-channel  is  equal  and, 
once  transparent,  no  further  changes  to  the  relative  phase 
or  amplitude  of  each  field  occur  [7].  Additionally,  we  ob¬ 
serve  important  saturation  mechanisms  which  may  limit 
the  generated  radiation  intensity  attainable  with  this  type 
of  frequency  converter. 

We  begin  with  a  brief  theoretical  discussion  of  this  work. 
We  consider  a  four-state  system  in  the  rotating  wave  ap¬ 
proximation,  as  shown  in  the  inset  in  Fig.  1.  The  coupling 
laser,  tuned  to  the  |2)  — ►  |3)  transition  with  Rabi  frequency 
Oc  (406  nm),  and  the  probe  laser,  tuned  to  the  |1)  — * ►  |3) 
transition  with  Rabi  frequency  (283  nm),  adiabati- 
cally  establish  a  large  atomic  coherence  of  the  Raman 
transition,  which  is  described  by  the  off-diagonal  den¬ 
sity  matrix  element  p\2  [2].  A  third  laser  field  with  ra¬ 


dian  frequency  aje  is  tuned  to  resonance  with  the  |2)  — > 
|4)  transition  and  mixes  with  the  coherence  to  generate  a 
fourth  field  at  ojh  resonant  with  the  |1)  — +  |4)  transition. 
The  (small-signal)  electric-field  attenuation  coefficients  of 
the  mixing  and  generated  fields  are  ae  =  Naepnj2  and 
cth  =  No-hPn/2 ,  where  N  is  the  atom  density,  pjj  is 
the  total  fraction  of  atoms  in  state  |;),  and  cre  and  <x/, 
are  the  single-atom,  Doppler-broadened,  resonant,  absorp¬ 
tion  cross  sections.  The  coherently  trapped  population 
is  distributed  according  to  pi\  =  SI* /(SI*  +  SI*),  pn  — 
ftJ/( n2c  +  ftp),  andpl2  =  -ftpft*/(ft2  +  Q2)  With 
these  definitions,  in  a  plane-wave  approximation  and  ignor¬ 
ing  saturation  effects  and  light  shifts,  the  propagation  equa¬ 
tions  for  a  quadrupiy  resonant  frequency  converter  from 
(le  (233  nm)  to  Sih  (186  nm)  are  [2] 


FIG.  1.  Schematic  diagram  of  the  experimental  apparatus. 
BC1,  etc.,  described  in  main  text.  Inset:  Energy-level  diagram 
of  the  double-A  system  in  208 Pb  vapor.  Laser  fields  with 
Rabi  frequencies  Slp ,  Oc,  and  Sle  are  applied,  and  Slh  is 
generated.  State  designations  and  energies:  |1),  6 p2  3/>0, 
0.0  cm"1  (ground);  |2),  6 p2  2Pj,  10650.3271  cm'1;  |3),  6pls 
3P,,  35  287.2244  cm"1;  |4>,  6p9s  3P,,  53 51 1.1485  cm”1.  The 
ionization  potential  of  Pb  I  is  59  819  cm”1  (7.42  eV). 
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+  aene  =  p\2Clht\ p(+i&kz),  (la) 


—  +  ahtlh  =  pu n,exp(-»A*z).  (lb) 

dz  2 

Since  ^-vector  matching  of  the  interacting  laser  fields  is 
automatically  satisfied  by  resonant  copropagating  fields  in 
the  rotating  wave  approximation,  Ak  =  0.  The  solutions 
of  Eqs.  (1)  with  boundary  conditions  fl/,(z  =  0)  =  0  and 
(1*0:  =  0)  —  fle(0)  are 


nh(z) 

Cle(  0) 


iNcrhp\2  exp 


ajz  ]  sin(/arz/2) 

2  J  aj 


£ ]eiz) 
0,(0) 


(£) 

(-) 

(1  -  e~a'z). 

(2a) 

v  nc  / 

\ar  / 

— —  + 

ae 

—  e 

~aTz 

(2b) 

aT 

OCj 

where  aT  =  ote  +  a*.  This  solution  assumes  that  the 
atomic  coherence  p\i  established  by  the  fields  in  the  first 
A-channel  is  unaffected  by  the  generation  process;  this  is 
ensured  by  maintaining  0,(0)  an  order  of  magnitude  less 
than  Op  and  O,  during  the  experiment. 

The  central  result  of  this  work  is  that,  in  the  limit 
of  many  absorption  lengths  ( aTz  — *  x),  the  nonlinear 
generation  process  matches  the  ratio  of  Rabi  frequencies 
in  each  A-channel  so  that  0/,(z)/0,(z)  — *  Op/Oc.  The 
frequency  conversion  behavior  predicted  by  Eq.  (2)  is 
dramatically  different  from  that  witnessed  in  earlier,  far- 
detuned  experiments  [1];  no  longer  is  there  parametric 
power  oscillation  between  fields  as  a  function  of  distance 
into  the  nonlinear  medium,  but,  rather,  monotonic  growth 
towards  an  asymptotic  value.  This  result  may  be  viewed 
as  a  special  case  of  matched  pulses  [7]  in  a  double-A  sys¬ 
tem;  atoms  evolved  into  the  coherent  population-trapped 
eigenstate  are  invisible  to  any  pair  of  optical  frequency 
fields  on  the  |1)  |4>  and  |2)  |4)  transitions  whose 

complex  ratio  of  Rabi  frequencies  is  flp/ilc.  Pairs  of 
input  fields  which  differ  from  this  ratio  undergo  selec¬ 
tive  absorption  of  unmatched  Fourier  components  and 
nonlinear  generation  until  they  conform  to  this  ratio  and 
are  thereby  decoupled  from  the  medium.  This  behavior 
is  quite  different  from  the  density  invariance  of  previous 
resonant  frequency  conversion  experiments  caused  by 
competition  between  four-wave  mixing  and  two-photon 
absorption  [8]. 

We  now  turn  to  a  description  of  the  experiment.  A 
schematic  diagram  of  the  experimental  setup  is  shown  in 
Fig.  1.  Three  linearly  polarized,  coherent  beams  at  406, 
283,  and  233  nm  are  spatially  and  temporally  overlapped 
by  beam  combiners  BC1  and  BC2,  converted  to  circu¬ 
lar  polarization  by  a  Fresnel  rhomb  FR,  and  propagate 
collinearly  through  the  Pb  atomic  vapor  contained  in  a 
heated  sidearm  cell.  The  output  beams  propagate  ~-T  cm 
in  atmosphere  into  a  vacuum  chamber  VC  (maintained 
<500  mTorr)  and  are  dispersed  by  an  isosceles  Brewster 
prism  IBP.  The  generated  field  at  186  nm  and  residual 
mixing  field  at  233  nm  are  directed  onto  sodium  salicy¬ 


late  scintillators  S.  The  fluorescence  signals  are  sampled 
by  fast  photodiodes  PD  l  and  PD2.  Each  waveform  is 
digitized  and  downloaded  to  a  computer.  The  data  are 
analyzed,  sorted  by  the  pulse  timing,  and  plotted;  data  with 
pulse  timing  overlap  errors  less  than  5  ns  are  retained.  The 
data  shown  in  Figs.  2-5  are  from  individual  events  and  are 
not  averaged. 

The  laser  setup  is  similar  to  that  used  in  previous 
experiments  [2,9].  All  beams  are  collimated  with  in¬ 
tensity  FWHM  beam  diameters  of  2.0  mm  (coupling 
laser,  406  nm,  30-ns  pulse  width),  2.12  mm  (probe  laser, 
283  nm,  15-ns  pulse  width),  and  1.24  mm  (mixing  laser, 
233  nm,  15-ns  pulse  width).  The  area  of  the  generated 
beam  is  inferred  to  be  the  same  as  that  of  the  mixing  laser. 
Laser  frequencies  are  determined  using  a  Burleigh  Model 
4500  pulsed  wave  meter  with  a  maximum  resolution  of 
0.01  cm'1  (300  MHz).  Typical  energies  and  (peak)  Rabi 
frequencies  for  the  probe,  coupling,  mixing,  and  generated 
lasers  are  50  /xJ  (0.25  cm-1),  50  pi  (0.30  cm-1),  1  pi 
(0.0148  cm-1),  and  0.1  pi  (5.31  X  10~3  cm-1),  respec¬ 
tively.  (These  beam  energies  are  a  factor  of  200  lower 
than  those  employed  in  previous  EIT-enhanced  nonlinear 
optics  experiments  [1,2].) 

The  186-nm  vacuum  ultraviolet  radiation  is  generated  as 
a  collimated  beam  which  propagates  collinearly  (i.e.,  on 
axis)  with  the  three  applied  fields.  Typical  pulse  energies 
are  less  than  1  pi  and  are  determined  using  a  scintillator/ 
photodiode  combination  which  has  been  calibrated  with 
known  233-nm  beam  energies  [10]. 

The  atomic  vapor  is  99.86%  isotopically  enriched  208 Pb 
metal  (Oak  Ridge  National  Laboratory),  heated  within  a 
25-cm-long  fused  quartz  sidearm  cell  manufactured  by 
Opthos  Instruments.  As  in  previous  work  [11],  densities 


FIG.  2.  Ratio  of  generated  186-nm  and  output  233-nm  Rabi 
frequencies,  at  the  exit  of  the  cell  (1*0  =  L)/Q,e(z  =  L ),  as 
a  function  of  the  coupling  laser  Rabi  frequency,  at  fixed 
(a)  For  NL  product  of  1.2  X  10 14  cm-2,  (b)  For  NL  product  of 
5.6  X  1014  cm-2.  The  dashed  vertical  line  in  each  figure  is  the 
value  of  Hc  required  for  weak-probe  EIT.  The  hyperbola  is  the 
aTL  — ►  theoretical  prediction  of  Eq.  (2),  taking  equal  to 

(a)  0.25  cm-1  and  (b)  0.24  cm-1. 
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FIG.  3.  Ratio  of  generated  186-nm  and  output  233-nm  Rabi 
frequencies  at  the  exit  of  the  cell  Qh(z  —  L)/(le(z  =  Lk  as 
a  function  of  the  probe  laser  Rabi  frequency,  at  fixed  f lc. 

(a)  For  NL  product  of  1.2  X  1014  cm-2,  (b)  For  NL  product 
of  5.6  X  1014  cm  2.  The  solid  line  is  the  <xjL  — *  sc  theoretical 
prediction,  taking  equal  to  (a)  0.33  cm-1  and  (b)  0.35  cm"1. 

are  determined  using  the  curve-of-growth  method  (to  an 
accuracy  of  ±10%).  At  600  °C,  the  resonant  absorp¬ 
tion  cross  sections  are  cre  =  2.21  X  10'13  cm2  and  cr^  — 
2.83  X  10  13  cm2.  We  investigate  the  behavior  of  the 
quadruply  resonant,  frequency  converter  at  two  density- 
length  (NL)  products:  (a)  1.2  X  1014  cm’2  (588  °C)  and 

(b)  5.6  X  1014  cm'2  (650  °C). 

In  Fig.  2,  we  study  the  dependence  of  the  ratio  of 
the  Rabi  frequencies  of  the  generated  VUV  and  output 
mixing  fields  at  the  exit  of  the  sidearm  cell,  =  L)/ 
Cle(z  =  L),  on  the  coupling  laser  Rabi  frequency  Ilc. 
Here,  the  probe  laser  Rabi  frequency  is  held  at  (a)  0.25  ± 


FIG.  4.  Ratio  of  the  output  and  input  233-nm  Rabi  frequencies 
fL(c  =  L)/de(z  =  0)  as  a  function  of  the  ratio  of  probe  and 
coupling  laser  Rabi  frequencies  dp/dc.  (a)  For  NL  product  of 
1.2  X  1014  cm-2,  (b)  For  NL  product  of  5.6  X  1014  cm"2.  The 
solid  line  is  the  aTL  — ►  x  theoretical  prediction. 


0.01  cm'1  and  (b)  0.24  ±  0.01  cm"1,  and  He(0)  < 
0.025  cm'1.  The  curved,  solid  lines  represent  the 
theoretical  prediction  of  Eqs.  (2)  at  these  values  of 
ap,  assuming  complete  Rabi-frequency  matching  (i.e., 
taking  arL  —*  “).  The  dashed,  vertical  line  indicates 
the  coupling  laser  Rabi  frequency  at  which,  for  each 
density,  the  step-up  in  EIT  occurs  (i.e.,  where,  in  the 
absence  of  the  mixing  laser,  the  transmission  of  a  weak 
probe  laser  field  through  the  cell  exceeds  50%).  Above 
this  value,  sufficient  preparation  energy  has  been  applied 
to  the  system  to  ensure  evolution  of  the  ensemble  into 
the  population-trapped  eigenstate  [12].  In  the  following 
figures,  only  those  points  with  Clc  >  0.3  cm'1  are  shown. 

In  Fig.  3  we  show  (lh(z  =  L)/de(z  =  L)  as  a  function 
of  the  probe  laser  Rabi  frequency  dp  at  fixed  dc  and  the 
ajL  — >  theoretical  prediction.  A  striking  feature  of 
these  data  is  the  departure  from  theory  when  the  Rabi  fre¬ 
quency  of  the  probe  laser  exceeds  that  of  the  coupling  laser. 
This  is  most  likely  related  to  a  failure  of  EIT;  increasing 
the  probe  laser  strength  relative  to  the  coupling  laser  in  the 
presence  of  nonzero  Raman-transition  dephasing  rates  ex¬ 
acerbates  two-photon  absorption  and  causes  non-negligible 
saturation  and  loss  of  coherence  on  time  scales  com¬ 
parable  to  the  pulse  widths.  (The  dephasing  time  T2  of 
the  Raman  transition  coherence  is  ^30  ns.)  We  expect 
that  the  use  of  shorter  (but  still  adiabatic)  pulses  in  such 
an  experiment  would  limit  this  saturation  and  increase 
the  strong-probe  conversion  efficiency.  This  strong-probe 
saturation  behavior  has  been  observed  in  other  work  in  a 
continuous-wave  regime  [13].  We  limit  display  of  data  in 
the  remaining  figures  to  those  points  where  Clp  <  dc. 

In  Fig.  4,  we  study  the  ratio  of  output  and  input  mixing 
field  Rabi  frequencies  fMz  =  L)/f!f(0)  as  a  function  of 


FIG.  5.  Evolution  to  matched  pulses.  Ratio  of  generated 
186-nm  and  output  233-nm  Rabi  frequencies  at  the  exit  of 
the  cell  dh{z  =  L)/de(z  =  L),  as  a  function  of  the  ratio  of 
probe  and  coupling  laser  Rabi  frequencies  dp/Slc.  (a)  For 
NL  product  of  1.2  X  1014  cm"2,  (b)  For  NL  product  of 
5.6  X  1014  cm"2.  The  solid  line  is  the  aTL  — ►  »  theoretical 
prediction. 
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the  ratio  of  the  probe  and  coupling  laser  Rabi  frequencies 
and  compare  with  the  predictions  of  Eq.  (2b).  Variation 
with  the  probe  and  coupling  lasers  is  caused  by  the  im¬ 
plicit  dependence  of  the  attenuation  coefficients  ae  and 
a:/,  on  the  population  terms  pn  and  p\\.  Note  that,  in  the 
limit  of  weak  probe  (vanishing  coherence),  no  population 
is  transferred  to  state  |2),  and  there  is  no  appreciable  ab¬ 
sorption  of  the  mixing  field. 

The  data  of  the  previous  figures  are  combined  in  Fig.  5 
to  illustrate  the  evolution  to  a  Rabi-frequency  matching 
regime  at  the  higher  temperature.  The  opacity  of  the 
low  density  case  (arL  ~  13)  should  be  more  than  suf¬ 
ficient  to  establish  the  Rabi-frequency  matching  condi¬ 
tion,  and  no  density  dependence  should  remain.  It  is  clear 
from  the  data,  however,  that  matching  is  not  attained  un¬ 
til  the  higher  density  case  (arL  =  57).  The  discrepancy 
arises  from  neglecting  saturation  and  power  broadening 
of  the  \l)  *-+  |4)  and  |2)  <-►  |4)  transitions  in  the  deriva¬ 
tion  of  Eqs.  (1).  This  assumption  requires  that  the  mix¬ 
ing  and  generated  field  Rabi  frequencies  f le  and  fi/,  are 
both  much  less  than  the  natural  linewidth  of  the  transi¬ 
tions.  During  the  experiment,  the  intensities  of  these  fields 
are  usually  an  order  of  magnitude  larger  than  the  satu¬ 
ration  intensities.  We  identify  power  broadening  by  the 
233-nm  field  and  the  resulting  decrease  in  the  number  of 
absorption  depths  at  both  oje  and  u)h  as  the  primary  satu¬ 
ration  mechanism  in  this  experiment.  Close  agreement 
between  numerical  simulation  and  experimental  results  is 
obtained  if  the  small-signal  absorption  coefficients  ae  and 
an  are  replaced  by  their  power-broadened  values.  Satura¬ 
tion  is  most  evident  in  the  behavior  of  the  conversion  effi¬ 
ciency  from  (De  — • ►  u)h  \  small-signal  efficiencies  [fMO)  < 
0.005  cm-1]  exceed  30%  but  drop  sharply  to  between 
5%  and  10%  when  fi^(0)  is  raised  to  0.02  cm-1  because 
the  effective  opacity  is  progressively  reduced  to  the  point 
where  the  matching  condition  is  no  longer  attained. 

At  the  higher  NL  condition,  we  achieved  approximately 
1%  overall  energy  conversion  efficiency.  A  VUV  pulse 
energy  of  1  p.J  was  generated  using  10  /xJ  of  233  nm  and 
50  /xJ  each  of  283  and  406  nm. 

In  summary,  we  have  demonstrated  efficient  nonlinear 
frequency  conversion  in  a  quadruply  resonant  double-A 
atomic  system,  using  collimated  laser  fields  with  modest 
power  densities.  This  work  provides  the  first  experimen¬ 
tal  realization  of  nonlinear  optics  in  a  density-invariant, 
Rabi-frequency  matching  regime.  Three  resonant,  arbi¬ 
trarily  phased  fields  are  applied  at  z  ~  0,  and  after  a  char¬ 
acteristic  distance,  the  atom  renders  itself  transparent  to 
all  fields  through  generation  of  a  properly  phased,  fourth 
resonant  field.  Transparency  is  achieved  once  the  com¬ 
plex  ratio  of  Rabi  frequencies  in  each  A-channel  is  equal, 
and  once  transparent,  no  further  changes  to  relative  phase 
or  amplitude  occur.  The  overall  conversion  efficiencies  of 
these  all-resonant  frequency  converters  are  limited  by  the 
preparation  energy  requirement  necessary  to  establish  EIT 
in  the  first  A-channel  and  by  power  broadening  of  the  sec¬ 
ond  A-channel.  We  note  that  the  use  of  a  single-isotope 


atomic  medium  minimizes  the  required  values  of  Clc  [14] 
but  is  otherwise  unnecessary. 
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Note  added  in  proof — A  comprehensive  theoretical 
study  of  nonlinear  optics  in  double-A  atomic  systems, 
with  emphasis  on  pulse-matching  properties,  has  recently 
been  published  by  Korsunsky  et  aL  [15]. 
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We  use  molecular  deuterium  and  two  driving  lasers  to  demonstrate  collinear  generation  of  mutually 
coherent  equidistant  sidebands,  covering  50000  cm-1  of  spectral  bandwidth  and  ranging  from  2.94  /xm 
to  195  nm  in  wavelength.  The  essential  idea  is  the  adiabatic  preparation  of  a  single,  highly  coherent 
(\pab\  —  0.33)  molecular  eigenstate. 


PACS  numbers:  42.50.Gy,  32.80.Qk,  42.60.Fc,  42. 65. Dr 

It  has  recently  been  predicted  that  coherent  molecular 
oscillations  can  produce  laser  frequency  modulation  (FM), 
with  a  total  bandwidth  extending  over  the  infrared,  visible, 
and  ultraviolet  spectral  regions,  and  with  a  possibility  of 
subfemtosecond  pulse  compression.  The  essence  of  this 
technique  is  the  use  of  a  Raman  transition  with  a  suffi¬ 
ciently  large  coherence  that  the  generation  length  and  the 
phase-slip  length  are  of  the  same  order.  This  coherence 
is  established  by  driving  the  molecular  transition  with  two 
single-mode  laser  fields,  slightly  detuned  from  the  Raman 
resonance  so  as  to  excite  a  single  molecular  eigenstate. 
Molecular  motion,  either  in  phase  with  the  driving  force 
(Raman  detuning  below  resonance)  or  antiphased  (detun¬ 
ing  above  resonance),  in  turn  modulates  the  driving  laser 
frequencies,  causing  the  collinear  generation  of  a  very 
broad  FM-like  spectrum  [1]. 

This  Letter  describes  the  first  experimental  results  for 
this  technique.  Using  molecular  deuterium  (£>2),  we 
demonstrate  collinear  generation  of  a  Raman  spectrum 
extending  over  50000  cm-1.  The  spectrum  consists  of 
seventeen  mutually  coherent  sidebands,  spaced  by  the 
fundamental  vibrational  frequency  of  £>2  (2994  cm-1), 
and  ranging  in  wavelength  from  2.94  /xm  to  195  nm. 

There  is  prior  work  on  broadband  coherent  Raman 
generation.  Several  groups  have  studied  the  enhancement 
of  stimulated  Raman  scattering  by  application  of  the  first 
Stokes  component  collinearly  with  the  pump  laser  beam; 
they  have  demonstrated  increased  conversion  efficiency, 
improved  reproducibility,  and  small  divergence  of  the 
anti-Stokes  beams  [2-4].  Hakuta  etal  have  demonstrated 
collinear  Raman  sideband  generation  in  solid  molecular 
hydrogen  (PL)  [5].  Nazarkin  et  al  have  shown  efficient 
laser  frequency  modulation  by  impulsively  exciting  a 
coherent  vibration  of  SF6  [6].  In  other  experiments, 
Ruhman  et  al.  have  used  impulsive  stimulated  Raman 
scattering  to  observe  coherent  molecular  vibrations  in  the 
time  domain  [7],  The  phased  and  antiphased  excitation 
techniques,  and  the  resultant  maximum  coherence  of  a 
molecular  state,  have  not  been  discussed. 

Our  experimental  setup  is  shown  in  Fig.  1.  We  use 
two  transform-limited  laser  pulses  at  wavelengths  of 
1.0645  /xm  and  807.22  nm,  such  that  the  (tunable)  laser- 
frequency  difference  is  approximately  equal  to  the  funda¬ 
mental  vibrational  frequency  in  £)2.  The  first  laser  is 

0031-9007/ 00/85(3)/562(4)$15.00 


a  Quanta-Ray  GCR-290  Q-switched  injection-seeded 
Nd:YAG  laser.  Its  output  is  attenuated  to  produce  100  mJ, 
12  ns  transform-limited  pulses  at  a  10  Hz  repetition  rate. 
The  laser  linewidth  is  calculated  to  be  <5  wiaser  =  37  MHz. 
The  second  laser  is  a  lab-built  Ti:sapphire  system,  injection 
seeded  from  an  external-cavity  laser  diode  and  pumped  by 
the  second  harmonic  of  a  separate  Q-switched  Nd:YAG 
laser.  This  laser  produces  75  mJ,  16  ns  transform-limited 
pulses  at  the  seeding  laser  wavelength.  This  wavelength 
can  be  tuned  precisely  and  is  monitored  by  a  Burleigh 
WA-1500  Wavemeter  (resolution  of  50  MHz).  The  two 
driving  laser  pulses  are  synchronized  by  adjusting  the 
delay  between  the  two  Nd:YAG  laser  Q-switch  trigger 
pulses.  The  laser  beams  are  combined  on  a  dichroic  beam 
splitter  and  are  loosely  focused  to  a  nearly  diffraction-lim¬ 
ited  spot  in  a  D2  cell.  The  1.06  /xm  laser  spot  size  is 
460  /xm,  and  the  807  nm  laser  spot  size  is  395  /xm. 

We  determine  the  Raman  resonance  by  scanning  the 
Ti:sapphire  laser  frequency  and  measuring  weak  Stokes 
gain  at  a  low  pump  intensity  and  a  small  D2  density. 
The  D2  cell  is  cooled  by  liquid  nitrogen  to  T  =  77  K; 
the  length  of  the  cooled  region  is  50  cm.  Cooling  re¬ 
duces  the  Doppler  linewidth  to  260  MHz  and  increases 
the  population  of  the  ground  rotational  state  of  D2  to  60%. 


FIG.  1.  Experimental  setup  and  energy  level  diagram  for  co¬ 
herent  molecular  excitation  and  collinear  Raman  generation.  The 
Raman  detuning  Aw  (positive  as  shown)  is  set  by  the  driving 
laser  frequencies. 
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The  collisional  dephasing  iinewidth  in  D 2  at  T  =  77  K  is 
450  kHz/torr.  as  deduced  from  room-temperature  data  [8], 

When  we  tune  the  driving  infrared  lasers  to  within 
1  GHz  from  the  Raman  resonance,  we  see  a  bright  beam 
of  white  light  at  the  output  of  the  D-_  cell.  We  disperse  the 
spectrum  with  a  prism  and  project  (without  collimating 
or  focusing)  onto  a  white  scintillating  screen  (60  cm 
wide,  270  cm  away  from  the  cell).  At  a  £>2  pressure  of 
50-100  torr  and  at  a  Raman  detuning  of  a  few  hundred 
MHz,  we  observe  up  to  thirteen  anti -Stokes  sidebands 
and  two  Stokes  sidebands,  in  addition  to  the  two  driving 
frequencies.  These  sidebands  are  spaced  by  2994  cm-1 
and  range  from  2.94  gm  to  195  nm  in  wavelength.  The 
next  Stokes  sideband  has  a  wavelength  of  24  /xm  and  is 
difficult  to  detect.  The  next  (fourteenth)  anti-Stokes  side¬ 
band  has  a  wavelength  of  184  nm  and  is  absorbed  by  air. 

Figure  2  shows  pictures  of  the  spectrum  taken  with 
a  digital  color  camera  with  a  single-shot  exposure  at  a 
fixed  aperture  size.  Starting  from  the  left,  the  first  two 
sidebands  are  the  driving  frequencies,  and  the  next  four 
are  anti-Stokes  sidebands  in  real  color  (red,  green,  blue, 
and  violet);  beginning  at  the  fifth  anti-Stokes,  the  side¬ 
bands  are  in  the  ultraviolet  and  only  fluorescence  is  visible. 
Figures  2(a)— 2(c)  show  the  spectrum  generated  at  a  £>2 
pressure  of  P  —  71  torr  and  a  Raman  detuning  of  Acu  = 
-400  ±  25  MHz  in  . part  (a),  Ao>  =  100  ±  25  MHz  in 
part  (b),  and  Acu  =  700  ±  25  MHz  in  part  (c).  The 
smooth  near-Gaussian  beam  profiles  for  nearly  all  side¬ 
bands,  as  shown  in  Figs.  2(a)— 2(c),  demonstrate  collinear 
anti-Stokes  generation  in  a  regime  of  high  molecular  co¬ 
herence.  At  higher  pressures  the  generation  is  no  longer 
collinear  and  the  anti-Stokes  sidebands  emerge  in  circles  of 
increasing  diameter.  An  example  at  a  pressure  of  350  torr 
and  Aa>  =  700  MHz  is  shown  in  Fig.  2(d). 

Figure  3  shows  Stokes  and  anti-Stokes  energy  spec¬ 
tra  for  three  different  values  of  the  Raman  detuning  and 
contrasts  the  on-resonance  generation  at  P  =  72  torr  (tri¬ 
angles  and  solid  line)  with  the  generation  below  resonance 


(circles  and  dashed  line)  and  above  resonance  (squares 
and  dotted  line).  We  observe  that  on-resonance  genera¬ 
tion  is  less  efficient  than  off-resonance  generation  for  all  of 
the  anti-Stokes  sidebands,  and  that  generation  below  reso¬ 
nance  is  more  efficient  than  above  resonance.  Note  that  we 
generate  hundreds  of  yuJ  per  pulse  (at  a  10  Hz  repetition 
rate)  at  sidebands  far  into  the  UV. 

We  now  proceed  with  the  discussion  of  these  results. 
We  consider  a  set  of  equidistant  Raman  sidebands  with 
complex  field  envelopes  Eq  and  carrier  frequencies  toq , 
tuned  close  to  the  Raman  transition  | a)  — *  \b ),  and  far 
detuned  from  the  upper  electronic  states  |i).  Propagation  of 
each  sideband  is  governed  by  the  slowly  varying  envelope 
equation  in  local  time  [1]: 


j  T)H(OqN  {ClqPaaEq  +  dqp^Eq  4-  b  ^pa^Eq~\ 


+  b<t+\PabEq+0*  0) 

where  paa  and  p^  are  populations  of  states  |a)  and  | b). 
Pub  is  the  Raman  coherence,  N  is  the  molecular  density, 
and  7]  =  (/jl/sq)1''2.  Expressions  for  the  dispersion  and 
coupling  constants  aq ,  bq,  and  dq  are  given  in  Ref.  [1], 

When  the  magnitude  of  the  coherence  pa\}  approaches 
its  maximum  value  of  0.5,  the  coupling  terms  in  Eq.  (1 ) 
(last  two  terms)  are  of  the  same  order  as  the  dispersion 
terms  (first  two  terms).  In  such  a  high-coherence  regime, 
the  generation  length  becomes  comparable  to  the  phase- 
slip  length,  and  efficient  generation  proceeds  collinearly, 
with  phase  matching  playing  a  small  role  [9].  To  the  ex¬ 
tent  that  dispersion  can  be  neglected  completely,  theory 
predicts  a  pure  sinusoidal  frequency  modulation  of  a  laser 
beam,  as  a  result  of  its  propagation  through  the  coherent 
medium  [1], 

The  excitation  of  the  molecular  states  | a)  and  \b)  is 
described  by  an  effective  two-by-two  Hamiltonian  [1]: 


tfeff  =  - 


A 

B* 


B 

D  -  2Ao>  ’ 


(2) 


FIG.  2  (color).  Spectrum  generated  in  the  setup  of  Fig.  1  at  (a)  P  =  71  torr  and  Aa>  —  —400  MHz.  (b)  P  —  71  ton  and  A u>  — 
100  MHz,  (c)  P  =  71  ton  and  lo)  =*  700  MHz,  and  (d)  P  =  350  ton  and  A  id  =  700  MHz.  We  observe  the  two  driving  infrared 
fields  (on  the  left),  and  multiple  generated  visible  and  ultraviolet  anti-Stokes  sidebands.  To  reduce  camera  saturation  for  (a)-(c),  the 
first  four  anti-Stokes  beams  are  attenuated  by  a  factor  of  100;  the  driving  field  at  807  nm  and  also  the  fifth  anti-Stokes  beam  are 
attenuated  by  10.  In  (d),  the  807  nm  beam  and  the  first  four  anti-Stokes  beams  are  attenuated  by  10. 
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FIG.  3.  Pulse  energies  generated  in  the  setup  of  Fig.  1  at  P  — 
72  torn  The  triangles  show  on-resonance  generation  (Aa>  =  0), 
the  circles  show  generation  by  phased  (Aa>  *=  500  MHz)  and 
the  squares  by  anti-phased  (Aw  =  -200  MHz)  states  of  Di. 

where  A  =  Y,qaq\Eq\2i  B  —  Y*qbqEqEq-\,  and  D  — 
£qdq\Eq\2.  The  Raman  detuning  Aco  is  the  difference 
between  the  molecular  transition  frequency  (Oy  —  coa 
and  the  frequency  of  the  driving  force  coq  -  a>q-\.  The 
collisional  dephasing  is  included  in  the  equation  for  the 
off-diagonal  density-matrix  element  (the  coherence)  pab. 
For  the  preparation  of  maximum  coherence  it  is  essential 
that  the  dephasing  time  is  longer  than  the  laser  pulse 
length  (as  is  the  case  in  our  experiment). 

The  molecular  eigenstate  which  evolves  from  the  ground 
state  as  the  field  amplitudes  are  increased  and  the  coher¬ 
ence  which  is  obtained  from  it  are 

|+)  =  cosy  exp^y  y)l«>  +  siny  exp^-y  y-j|i>, 

!  (3) 

Pab  =  y  sin0exp(y», 

where  B  =  |S|exp(y‘(p)  and  tan#  =  2|f?|/(2Aru  —  D  + 
4).  When  the  Raman  detuning  is  larger  than  the  driving 
laser  linewidth,  the  molecular  system  follows  this  eigen¬ 
state  adiabatically,  with  the  sign  of  pab  determined  by  the 
sign  of  the  Raman  detuning.  In  a  resonant,  nonadiabatic 
regime  the  molecules  exhibit  two-photon  Rabi  flopping  at 
a  frequency  B. 

In  our  experiment  the  Doppler  broadened  linewidth 
(260  MHz)  of  the  2994  cm"1  transition  of  D2  is  larger 
than  both  the  collisional  dephasing  rate  (32  MHz  for 
P  =  72  torr)  and  the  laser  linewidth  (37  MHz).  When 
we  tune  close  to  the  center  of  the  Doppler  line,  we 
expect  adiabatic  preparation  for  most  of  the  molecules  in 
both  wings  of  the  velocity  distribution,  with  significant 
magnitude  of  their  excitation  when  the  coupling  parameter 
B  >  5 laser*  The  phases  of  these  excitations  are  positive 
or  negative  depending  on  the  sign  of  the  Raman  detuning, 
such  that  the  contributions  of  these  molecules  partially 
cancel  each  other,  resulting  in  a  smaller  integrated  value 


of  the  coherence  pab  at  the  center  frequency.  However, 
when  we  detune  to  either  side  of  the  Doppler  line,  we 
expect  adiabatic  preparation  of  nearly  all  of  the  molecules 
in  the  same  (phase  or  antiphased)  state  [Eqs.  (3)],  re¬ 
sulting  in  a  larger  total  coherence  and  more  efficient 
Raman  generation.  Similar  behavior  has  been  observed  in 
saturated  coherent  anti-Stokes  Raman  spectroscopy  [10]. 

Figure  4  illustrates  this  behavior.  Parts  (a)  and  (b)  show 
the  first  anti-Stokes  pulse  energy  as  a  function  of  the  Ra¬ 
man  detuning,  both  at  low  pressure  (P  =  6  torr).  Here,  the 
depletion  of  the  applied  fields  is  not  important,  and  the  gen¬ 
erated  energy  is  proportional  to  the  square  of  the  Doppler- 
averaged  coherence  [Eq.  (1)].  In  Fig.  4(a),  the  1.06  fim 
laser  pulse  energy  is  5.4  mJ  and  B  —  39  MHz  ~  5<y|aser; 
for  these  conditions  the  anti-Stokes  generation  peaks  on 
resonance.  For  Fig.  4(b)  we  apply  full  power  of  the  driv¬ 
ing  lasers,  so  that  B  =  168  MHz  »  <5w]aser,  and  observe 
that  the  anti-Stokes  generation  peaks  on  both  sides  of  the 
resonance,  with  a  dip  exactly  on  resonance.  The  peak  at 
positive  A<w  is  higher,  because,  as  the  laser  power  increases 
as  a  function  of  time,  this  peak  Stark-shifts  toward  reso¬ 
nance.  We  solve  the  full  density-matrix  equations  numeri¬ 
cally  for  the  experimental  driving  laser  pulses,  calculate 
the  (complex)  coherence  pab  at  the  pulse  peak  for  different 
values  of  Raman  detuning,  and  average  it  over  the  Doppler 
distribution.  The  solid  lines  in  Figs.  4(a)  and  4(b)  show 
the  square  of  the  Doppler-averaged  coherence  (no  scaling). 
We  calculate  that  the  maximal  value  of  the  coherence  that 
we  achieve  in  our  experiment  is  Kpat>)DopPierl  =  0.33. 

Figure  4(c)  shows  qualitatively  similar  behavior  at 
a  higher  Di  pressure  (P  =  72  torr).  The  filled  circles 
represent  the  sum  of  all  generated  Stokes  and  anti-Stokes 
sidebands.  The  open  circles  and  triangles  show  the  first 
and  second  anti-Stokes  sidebands.  All  measurements  for 
Figs.  3  and  4  are  done  with  a  Molectron  J4-09  pyroelectric 
detector.  Each  point  in  the  plots  is  an  average  over  eight 
laser  shots,  with  error  bars  showing  the  standard  deviation 
for  each  data  set. 

In  order  to  investigate  adiabatic  and  nonadiabatic 
molecular  excitation  we  measure  the  total  laser  energy 
transmitted  through  the  Di  cell  (including  all  applied  and 
generated  sidebands)  as  a  function  of  the  Raman  detuning. 
When  the  molecular  evolution  is  adiabatic,  we  expect 
all  molecules  to  return  to  the  ground  state  as  the  laser 
fields  decrease  at  the  falling  edge  of  the  applied  pulses. 
We  observe  a  dip  in  transmission  on  Raman  resonance, 
confirming  that  energy  is  left  in  the  molecular  system  due 
to  nonadiabatic  behavior.  For  comparison  we  numerically 
calculate  the  probability  for  the  molecular  system  to 
remain  in  the  adiabatic  eigenstate  |+)  [Eqs.  (3)],  averaged 
over  the  Doppler  distribution  [solid  line  in  Fig.  4(d)], 
again  demonstrating  good  qualitative  agreement  between 
theory  and  experiment.  By  comparing  Figs.  4(c)  and  4(d), 
we  note  that  the  Raman  generation  peaks  at  detunings 
where  the  molecular  evolution  is  completely  adiabatic. 

We  note  the  connection  of  this  work  to  the  coherent 
population  transfer  [11]  and  to  electromagnetically  induced 
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FIG.  4.  Raman  generation  by  adiabatic  eigenstates  in  Do. 
(a)  First  anti-Stokes  generation  by  weakly  driven  vibrations 
( B  «=  laser)  at  P  =  6  torr.  (b)  First  anti-Stokes  generation 
by  strongly  driven  vibrations  (Z?  »  <5a>|aser)  at  P  =  6  torr.  In 
both  (a)  and  (b)  the  solid  lines  show  the  calculated  square  of 
the  Doppler-averaged  coherence  pab  at  the  pulse  peak,  (c)  Sum 
over  pulse  energies  of  all  generated  Raman  sidebands  (filled 
circles),  first  anti-Stokes  generation  (open  circles),  and  second 
anti-Stokes  generation  (triangles)  for  B  »  <5a>iaSer  and  P  = 
72  torr.  (d)  Total  energy  in  all  applied  and  generated  sidebands 
at  the  output  of  the  Di  cell.  The  solid  line  represents  the 
calculated  fraction  of  the  molecular  population  in  the  eigenstate 
which  evolves  adiabatically  from  the  ground  state  [Eq.  (3)]. 

transparency  (EIT)  [12]:  Each  involves  the  preparation 
of  the  adiabatically  evolving  eigenstate.  In  three-state  EIT 
this  is  done  by  using  quantum  interference  or  optical  pump¬ 
ing.  Here  it  is  done  by  fixed  or  adiabatically  changing 
detuning.  Both  EIT  and  this  process  require  that  the  num¬ 
ber  of  photons  in  the  laser  pulse  exceed  the  number  of 


molecules  in  the  laser  path  [12,13].  Both  allow  nonlinear 
optics  at  maximum  coherence,  where  the  role  of  phase 
matching  is  greatly  reduced  [9].  But,  differing  from  EIT, 
here,  the  refractive  indices  at  the  various  sidebands  are  re¬ 
duced  (in  the  antiphased  case)  [  l]  but  are  not  equal  to  unity. 

In  summary,  this  Letter  demonstrates  broadband 
collinear  Raman  generation  by  adiabatically  prepared 
phased  and  antiphased  states  in  molecular  deuterium. 
We  show  that,  in  agreement  with  theory,  generation 
maximizes  at  a  finite  detuning  on  either  side  of  the  Raman 
resonance.  We  believe  it  will  be  possible  to  recombine 
the  generated  sidebands  and  to  use  spectral  modification 
techniques  [14]  to  synthesize  specified  subfemtosecond 
time  structures  in  a  target  cell. 
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